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PREFACE 


Upon completion of my graduate work at the Naval 
Postgraduate School I will be assigned as an Assistant 
Professor in the Department of Mathematical Sciences, at the 
United States Military Academy. While there I will teach the 
cadets calculus. I will incorporate this package on integral 
calculus that I have developed into the cadets' instruction. 

As an instructor I will see this integral calculus package 
used by the cadets. I will get feedback from them on its use 
and its value in the instruction. I will then use the 
feedback to make improvements to the package and to expand its 
capabilities. It is my vision that such an instructional aid 
will improve the overall quality and effectiveness of the 
calculus instruction at the Military Academy. 
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ABSTRACT 


One of the most difficult concepts in the area of integral 
multivariate calculus is finding the limits of integration. 
This thesis describes an interactive computer program designed 
to help students understand this important concept. The 
program shows how a given domain is plotted and teaches how to 
find the limits of integration when evaluating two-dimensional 
integrals. 

The program allows the user to enter any known information 
about a region and then evaluates the integral. The region is 
plotted, the limits of integration are given along with the 
area of the region. The program handles cartesian and polar 
coordinate, two-dimensional integral problems. This program 
could be used independently by the student and/or used by the 
calculus instructor in the classroom. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 
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1. INTRODUCTION 


In the area of integral multivariate calculus one of the 
most difficult topics is deciding on the limits of integra¬ 
tion. This package shows the calculus student how a given 
domain is plotted and teaches the student how to find the 
limits of integration when evaluating two-dimensional 
integrals. 

Integral calculus is the mathematics we use, for example, 
to find lengths and areas of irregular shapes; to calculate 
average values of functions; and to predict future values in 
population sizes and future costs of living. The development 
of integral calculus starts from the calculation of areas. 
The limits used to define areas are special cases of a 
definite integral. 

The focus of the package is on teaching the student how to 
correctly find the limits of integration for a region. While 
for typical calculus problems the limits of integration are 
determined analytically, it is quite helpful to see the 
regions of interest. It then is important to be able to 
accurately graph the regions. The calculus student may have 
difficulty with curve sketching especially in polar 
coordinates. This package allows the student to visualize the 
region of interest enabling the student to better understand 
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the problem. It is essential that the student be able to see 
accurate plots. 

This package was developed using the Matbematica software. 
Matbematica uses an adaptive sampling algorithm for plotting 
to determine when and where to sample a function. A section 
with large curvature is sampled more frequently than a flat 
section of a function. This sampling algorithm coupled with 
a set prescribed number of equally spaced samples produces the 
quite reasonable plots we want. This graphing capability is 
one of the greatest attributes of this package. 

This package, which teaches the student how to find the 
limits of integration, has two parts. The first part, which 
is found in Appendix A, consists of a lab book. This lab book 
addresses two-dimensional integrals and walks the student 
through numerous example problems. The two-dimensional case 
considers regions given in the cartesian coordinate system as 
well as regions in the polar coordinate system. In all 
examples the student is shown the plot of the functions and 
then is shown how the limits of integration are determined. 
The integrals are set up and evaluated and the area of the 
region is given. 

The second part of this package is an interactive program 
written in Matbematica's programming language. With this 
program the student can solve two-dimensional integral 
problems, given in both the cartesian coordinate system and 
the polar coordinate system. The student enters the functions 
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which define a region's boundaries and any other known 
information about the region. The curves are then plotted. 
The limits of integration are determined and the area of the 
region is calculated. The program output includes the limits 
of integration, the curves and the calculated area of the 
regions. The program also can be used to display any of the 
plots for the example problems found in the lab book. 

The program's greatest value, however, is in allowing the 
user to input his own problems. The program is designed to 
handle those types of integral calculus problems found in any 
introductory level calculus text book. 

No knowledge of Mathematics is required to use the 
interactive program. The tutorial, found in Appendix B, tells 
the user how to load the program and how to use it. The 
tutorial takes the user through three example problems. The 
first is an integral problem given in cartesian coordinates. 
The second problem is in polar coordinates and finds the area 
between a polar curve and the origin. The final problem is 
also in polar coordinates and finds the area between two polar 
curves. These three problems are simple problems that show 
how the program is used. The tutorial, while not compre¬ 
hensive, gives the user some idea of the types of problems the 
program handles. It demonstrates how the program prompts the 
user for input and how that input is entered. The program is 
capable of much more complicated problems than those in the 
tutorial. In general, once the package is loaded and the 
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program is started the user simply follows the instructions 
that the program displ*.ys. 

Due to the variety of ways in which a region can be 
defined or in which a problem can be worded the program must 
be versatile. The program prompts the user in such a way as 
to Swiicit the necessary information about the region 
regardless of how the problem is stated. As mentioned 
previously the program is written to handle the types of 
problems found in calculus text books. These problems are 
generally well defined such that the program has no difficulty 
in determining the region of interest. On the other hand it 
is quite possible to contrive problems which are not well 
defined and for which the program is unable to determine the 
region. An example of this type is when there is not enough 
information provided. It is also possible to have a problem 
which is seemingly well defined for which the program does not 
work. An example of this type of problem is in cartesian 
coordinates when the functions intersect and do not define a 
single unique region. For instance the functions y=12-x^, 
y=x and y=0 intersect to form three separate regions (see 
Figure 1, page 5). 

If such a problem is entered, the program produces either 
no output, bad output, or an error message. There is still a 
way in which the program can be used to solve these problems. 
For problems given in cartesian coordinates there is a 
cartesian plot-only option. If this option is selected the 
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user can input all the functions for specified ranges of the 
variables. The functions will then be displayed and the plot 
can be used to estimate the limits of integration. With the 
estimated limits of integration obtained from the plot of the 
functions the user can run the program again with this new 
information. This plot-only option can be used also if the 
user just wants to see what a particular function looks like. 
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In general the program does not have any difficulty with 
the plots of the functions given in polar coordinates. If the 
limits of integration are not known in advance the program 
plots the polar functions from Q=0 to Q=2rr radians where Q 
represents the polar angle theta. For the two options, 
finding the area between a polar curve and the origin and 
finding the area between two polar curves, this feature acts 
as a plot-only option. The values of Q for which the polar 
functions intersect are displayed when finding the area 
between two curves. The values of Q for which the polar 
function is equal to zero are displayed when finding the area 
between a polar curve and the origin. 
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However if the user only w^nts to plot the polar functions 
without solving the integral problem, he can, both for 
plotting the area between a polar curve and the origin and for 
plotting the area between two polar curves. 
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II. MOTIVATION 


Integral calculus has a vide variety of applications in 
mathematics, physics, engineering, economics, medicine, and 
other fields. Integrals are used, for example, to compute 
areas, volumes, lengths of curves, forces of fluids, work, 
weight, centers of gravity, and probabilities, and to 
determine functions from their rates of change. In these many 
applications often the most difficult part is in properly 
setting up the integrals to be evaluated. More exactly the 
difficulty lies in determining the limits of integration for 
these integrals. 

Let R denote the region in the xy-plane that is bounded on 
the top by the graph y*f(x) and on the bottom by the graph 
y=g(x) and that extends from x=a on the left to x=b on the 
right. Double integrals over R may be expressed as iterated 
integrals by the formula 


( 1 ) 


x=b y=f{x) ^ 

JJ h(x,y)dA = J I J h(x,y)dy Jdx. 

R x=a L y=g(x) 

To evaluate the iterative integral we start at the innermost 

integral and work toward the outside. 

The most difficulty with evaluating two-dimensional 

integrals given in cartesian coordinates is in determining 

g(x) , f(x), a, and b which are limits of integration. The 

region R may have either changing upper or lower boundaries in 
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which case the region R is composed of subregions. The double 
integral over R then becomes equal to the summation of double 
integrals, that are over each of the subregions. Then for 
each subregion the limits of integration must be determined. 
This situation is more complicated, making determining the 
limits of integration difficult. 

The focus of this thesis is teaching how to find the 
limits. As such, all the two-dimensional integral problems in 
this thesis are to find the area of a given region. That is 
h(x,y)=l and so we have for integrals given in cartesian 
coordinates 


Area= 


x=b y=f(x) 

= J r J dy I dx. 

x=a y=g(x) 


( 2 ) 


After evaluating the innermost integral the equation becomes 


x=b 

Area= J [f(x)-g(x)Jdx. 
x=a 


(3) 


The functions y=f(x) and y=g(x) can be found by graphing 
the given functions. If the left limit, x=a, and the right 
limit, x=b, are not known they are determined analytically. 
The functions y=f(x) and y=g(x) are equated and solved for x. 
Similarly in polar coordinates we have 


rr r r^=^ 2 {Q) 

Area= JJ rdA = J J rdr dQ 

R Q=a r=rl(Q) 


(4) 
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where both rl and r2 are functions of Q. For finding the 
area between a polar curve and the origin, rl=0 and we have 


Q-b r=r2(Q) 

Area= J [ J rdr J dQ. 

Q*a 0 


If we evaluate the innermost integral we get 


.Q=b 2 

Area = J 1 r2 dQ. 
Q=a 2 


(5) 


( 6 ) 


For finding the area between two polar curves we evaluate 
the innermost integral in equation (4) above to get 

Q=b 2 2 

Area= J 1 [r2 - rl ] dQ. (7) 

Q=a ^ 

Again the most difficulty with two-dimensional integrals 
given in polar coordinates is setting up the integrals. More 
specifically the difficulty is in determining the elements a, 
b, rl, and r2 which are the limits of integration. In polar 
coordinate integral problems the graphs of the polar curves 
are used to find the limits of integration. In finding the 
area between a polar curve and the origin the limits a and b 
(if not known) are determined by solving the equation r2=0 for 
r. In finding the area between two polar curves the limits a 
and b (if not known) are determined by equating rl and r2 and 
solving for r. 





The calculus student must do numerous example problems in 
order to become proficient and confident in properly setting 
up an integral. The more problems a student sees and does, 
the more confident the student becomes. This package provides 
numerous worked solutions to problems. The program provides 
the solutions to as many integral problems as the student may 
wish to work. The program can be used with whatever calculus 
textbook the student may be using. 

The program gives solutions that include not only the 
numerical area, but also the graph of the region and the 
limits of integration. The students then can use this 
information to identify their problem areas. The students may 
find that they are getting the right limits of integration 
when they work the problems but that their solutions are not 
correct. Then it is likely that they are making mistakes when 
evaluating the integral. In this case the students see that 
they should review their integration techniques. They may be 
finding incorrect limits of integration. Then the graph of 
the region may be helpful in determining where the errors are 
being made. At any rate the program provides a great deal of 
useful information. Most textbooks provide just the value of 
the integral. 

The lab book contains numerous worked out example 
problems. It will complement the few worked problems that are 
likely to be found in most textbooks. The fact that the 
package is adaptable to any calculus textbook is very useful. 
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The simple-to-use program gives the students an extremely 
efficient and effective way of learning how to find limits of 
integration. The use of mathematical software can also make 
learning calculus more enjoyable. Typically students are 
first introduced to mathematical software in a higher level 
course. By using this package in their calculus instruction 
the students are introduced to mathematical software early 
during their education. In this age of computers this becomes 
more and more important all the time. 

At the Naval Postgraduate School there is currently no 
introduction to mathematical software included in the syllabus 
for the instruction of calculus. While it is true that the 
textbook being used references numerous examples and applica¬ 
tions of different mathematical software, typically no use of 
mathematical software is included in the course instruction. 

By introducing mathematical software in the instruction of 
calculus the students can be exposed to some of the very 
powerful uses and applications of mathematical software. With 
the students' introduction to mathematical software in their 
calculus instruction via this package a major step is being 
taken in changing current methods of instruction. This 
package could be used at the United States Military Academy as 
well as other colleges and universities. 
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III. USB OF PACKAOB IN CALCULUS INSTRUCTION 


This section addresses incorporating the integral package 
into the calculus instruction on two-dimensional integrals. 
The integral package is not intended to be used by itself to 
teach the calculus student how to find limits of integration 
and calculate two-dimensional integrals. Rather it should 
complement the classroom instruction and the homework problems 
on this topic. 

The package should be introduced to the students after 
they have received the classroom instruction on two- 
dimensional integrals given in cartesian coordinates. The 
students should also first attempt to work a number of 
assigned homework problems that require finding the area of a 
region given in cartesian coordinates. After the students see 
the solutions to these problems, they could be shown the 
package. 

The ideal way to introduce the package is in a computer 
lab in which all students can have hands-on use at a terminal. 
The instructor could either talk the students through turning 
on Mathematica, loading the package and using it or let the 
students simply follow the instructions in the tutorial found 
in Appendix B. At this point in the calculus instruction the 
students would only be ready for the first example problem in 
the tutorial. In this problem the students would use the 
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program to find the area of a region given in cartesian 
coordinates. 

The instructor could have other problems for the students 
to solve. These could be problems given in cartesian 
coordinates which come directly out of the calculus textbook 
being used. The greater variety of problems which the 
students solve using the program will demonstrate the 
program's many capabilities to solve integral problems given 
in cartesian coordinates. 

Once the students are shown how to use the integral 
program part of the package they should be introduced to the 
lab book which contains numerous worked problems. These 
problems show the students step-by-step how to find the limits 
of integration and evaluate the integral. These example 
problems will supplement the examples worked in the classroom 
instruction and those found in the calculus textbook being 
used. These detailed solutions will further reinforce the 
procedures established in the classroom for finding the limits 
of integration and setting up the proper integrals. 

The students at this point would not have received any 
classroom instruction on integrals involving polar 
coordinates. As such the students should only be told of the 
program's ability to handle such integral problems. The 
students should first be taught how to graph in polar 
coordinates. Once the students have received this classroom 
instruction on graphing and on integrals in polar coordinates 
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and had the opportunity to work some of these problems they 
should be shown the polar coordinate capabilities of the 
integral program. For polar coordinates the students should 
be shown how the program can find the area between a polar 
curve and the origin or to find the area between two polar 
curves. Again the tutorial is best suited for this, with one 
example of each of these problems. 

For all integral problems the student should first work 
the problems without the use of the program. The program 
should be used by the students to check their work and to 
identify any errors. All solutions given by the program 
include the plot of the functions and the region. The 
students can then see graphically presented the limits of 
integration. 

It is important that the student always first attempt to 
do the problems on their own. The program is not intended to 
do the student's homework but rather it is to assist the 
student in learning. The goal is that the students can, by 
themselves, find the limits of integration regardless how a 
problem is presented. The program is an educational aid and 
should be used as such. This integral package properly 
incorporated into the calculus instruction can make learning 
to find limits of integration for two-dimensional integrals 
easier and fun. This integral package complements the 
classroom instruction on finding the limits of integration for 
two-dimensional integrals. 
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IV. WHY MATHZMATICA 


vniile there are many types of mathematical software that 
I could have used, I selected Mathematiea. Since I have not 
used any other mathematical software extensively I am unable 
to comment on the relative advantages or disadvantages of one 
software over another. One of my primary reasons for 
selecting Mathematiea is its availability at both the Naval 
Postgraduate School and the United States Military Academy. 

Mathematiea also had a number of important attributes 
which lead me to select it. One such attribute of the 
software is its sophisticated graphing capabilities. 
Producing quality plots was essential to the package. 
Mathematiea enabled me to use a variety of graphing 
capabilities and options. Not only did they plot the 
functions but also provided necessary information about the 
regions. This information was used to find limits of 
integration, determine upper and lower curves, left and right 
curves, or inner and outer curves, and finally the area of the 
region. 

A second attribute of the software was its programming 
capability. Mathematics has an excellent pattern matching 
feature which was used extensively. I also found it easy to 
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program and debug in Mathematiea. With Professor David 
Canright's expertise in Mathematiea which was quite 
accessible, I felt confident in using it. 

Another advantage or rather attribute of Mathematiea is 
the ability for students to use the paclcage and other 
Mathematiea capabilities at the same time. In using 
Mathematiea I was also able to take advantage of many of 
Mathematiea'B built-in functions. 

Those platforms that support the program are IBM 
compatible personal computers with the Windows environment, 
the Unix workstation with X-Windows, and Macintosh computers. 
The program operation described herein is with the notebook 
interface and the Windows version of Mathematiea. The details 
differ for the non-notebook interface found with the Unix 
version of Mathematiea. 
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V. STRUCTfTRS OF THE PROGR1UC 


The integral prograa is organized into ten Nathefliatica 
packages which have nine major subprograms and numerous other 
defined functions. A package is simply a file containing 
Mathematlca definitions. To access the functions defined in 
the packages, the packages must first be loaded. The file 
named Integral.na is at the core of the integral program. 
This file connects the ten files together along with the 
functions defined within them. When integral.aa is loaded it 
in turns loads all the other files which make up the program. 

The program is started by calling the subprogram graph, 
found in file integral.ma. The rest of the program is built 
around subprogram graph which calls up the other subprograms. 
After the program is invoked subprogram graph prompts the user 
to respond to various instructions. It is in this manner in 
which different options are selected. Depending on which 
options are selected the appropriate subprograms are called 
(see Figure 2 on next page). 
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Program Architecture 


Subprogram 

graph 


To display a graph that 
is in the lab book type 
in the example number, 
otherwise type n. 


1 thru 23 



Type 1 for cartesian 
coordinate system, 
or type 2 for polar 
coordinate system. 


To just 
plot curves 
type y, 
otherwise 
type n. 




Type 1 for plot of area 
between the origin and a 
polar curve, or type 2 
for plot of area between 
two polar curves. 



Type 1 for the area 
between the origin and 
a polar curve, or type 2 
for the area between 
two polar curves. 






















If the user desires to see the plot of the region for any 
of the example problems in Appendix B, subprogram graph calls 
subprogram examples (found in file examples.me) . If the user 
just wants to plot a set of curves given in cartesian 
coordinates subprogram graph calls up subprogram cartplot 
(found in file cartplot.ma) . If the user has a cartesian 
coordinate problem to solve subprogram graph calls subprogram 
cart (found in file cart.ma) . The regions in these types of 
problems must be given as functions of the independent 
variable x. The integration is then with respect to x. For 
problems given in cartesian coordinates subprogram cart calls 
up subprogram Inverse (found in file invar so. ma) . Siibprogram 
inverse solves the problem by integrating over the region with 
respect to the variable y. 8\ibprogram inverse first 
determines if integration with respect to the variable y is 
appropriate. If the functions are invertible the user is 
then asked if he desires to integrate the same problem with 
respect to the variable y. 

If the user wants just a plot of the area between the 
origin and a polar curve, subprogram graph calls subprogram 
polplotl (found in file polplotl.ma) . If the user wants just 
a plot of the area between two polar curves, subprogram graph 
calls svibprogram polplot2 (found in file polplot2.ma) . 

If the problem is to find the area between a polar curve and 
the origin, subprogram graph calls subprogram polarl (found 
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in file polerl.M) . If the problem is to find the area 
between two polar cuzrves, subprogram graph calls subprogram 
polar2 (found in file polar2.ma). 

Many of the example problems in Appendix B have regions 
which are shaded. File shades.ma contains the Mathsmatlca 
definitions which produce the different shades for those 
example problems. 
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VI. SUBPROORAM GRAPH 


Subprograa graph is the main branch of the program. 
Depending upon the type of integral problem or option selected 
subprogram graph calls the appropriate subprograms. 

After the program is started the first prompt says: 

To display a graph that is in the lab book type in the 
example number, otherwise type n. 

The user can type a number from 1 to 23 to call up subprogram 
examples. Subprogram examples displays the graph which 
corresponds to that number example problem found in the lab 
book. 

If the user does not select this option the next prompt 
appears: 

Type 1 for cartesian coordinate system, or type 2 for 
polar coordinate systm. 

If the user types 1 for cartesian coordinate system the 
following prompt appears: 

To just plot curves type y, otherwise type n. 

If the user types y, subprogram graph calls up subprogram 
cartplot which is the cartesian plot-only option. Otherwise 
if the user types n, subprogram graph calls ups subprogram 
cart which solves those problems given in cartesian 
coordinates. 
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If when the user sees the prompt: 

Type 1 for oertesien coordinate system, or type 2 for 
polar coordinate system. 

He types 2 for polar coordinate system, the following prompt 
appears: 

To just plot curves type y, otherwise type n. 

If the user types y to select the polar plot-only option he 
then sees the next prompt: 

Type 1 for plot of area between the origin and a polar 
curve, or type 2 for plot of area between two polar 
curves. 

If the user types 1, subprogram graph calls up subprogram 
polplotl which allows the user to graph the area between a 
polar curve and the origin. Otherwise if the user types 2, 
sxibprogram graph calls up stibprogram polplot2 which allows the 
user to graph the area between two polar curves. 

If the polar plot-only option is not selected the user 
sees the following prompt: 

Type 1 for area between the origin and a polar curve, or 
type 2 for area between two polar curves. 

If the user types 1, subprogram graph calls up svtbprogram 
polar1 which finds the area bounded by a polar curve and the 
origin. Otherwise if the user types 2, subprogram graph calls 
up subprogram polar2 which finds the area between two polar 
curves. 
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VII. SUBPROGRAM CART 


Subprograa oart solves those problems given in cartesian 
coordinates. Subprogram graph calls up siibprogram cart when 
the user indicates he wants to solve a problem given in 
cartesian coordinates. Subprogram cart can be used by itself 
without invoking the main program. In order to use this 
stand-alone option the user types cart and then activates the 
cell. The tutorial covers the many ways in which the cell can 
be activated. 

Regardless how subprogram cart is called, the subprogram 
attempts to solicit all known information. The user is first 
told to enter the functions defining the region. These 
functions must be functions of the independent variable x. 
After the functions have been entered the user is told to 
indicate if any of the limits of integration is known. If the 
user indicates yes he is then prompted to input those known 
limits. The subprogram then includes the functions x=(left 
limit) and x=(right limit) as appropriate, to those functions 
of X already entered. If the left limit is not known and the 
right limit is either not known or greater than zero, the user 
is told to indicate if the region is restricted to the right 
half plane. If the right limit is not known and the left 
limit is either not known or less than zero, the user is told 
to indicate if the region is restricted to the left half 


23 





plane. In some problems the functions will intersect and form 
a region in the right half plane and one in the left half 
plane. The subprogram must know which region to consider. An 
example of this type of problem is the intersection of the 
functions y=12-x^, y=x and y=0 (see Figure 1 on page 5) . 

Once all this information is obtained, the subprogram is 
ready to solve the problem. It first obtains all the points 
of intersection for each of the curves and lines segments. 
For all pairs of equations it does this by simultaneously 
solving two equations with two unknowns, using Mathematica's 
NSolve function. The transcendental functions defined in the 
program are Taylor series expansions. When equations contain 
transcendental functions the subprogram uses these Taylor 
series expansions to find the approximate points of 
intersection. These approximations are then used as starting 
points in Mathematica's rindRoots function to find the points 
of intersection. FindRoots uses Newton's method to solve for 
these values. 

Subprogram cart uses the x coordinates of the points of 
intersection, for each function entered by the user, to define 
the right and left lijnits for the curve. If these x 
coordinates are the same value the subprogram knows that the 
points of intersection, for the entered function, form a 
vertical line segment. In order for the program to solve a 
problem each curve must have just two endpoints. 
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If there are fewer than two endpoints, the functions do 
not bound a closed region. If there are more than two points 
of intersection, the functions intersect to form more than one 
region. In the latter case, the subprogram is unable to solve 
the problem correctly. 

Consider the example where the lines y>x and y^sx are the 
only two functions entered. These functions intersect only at 
the point whose coordinates are (0,0). The program then 
determines only one endpoint (i.e., the point (0,0)) for each 
function. These two functions by themselves, do not bound a 
closed region and the program does not have enough informa¬ 
tion. 

Now consider the example where the functions are y=x, y=0 
and y=12-x^. The points of intersection, and thus the 
endpoints that the program determines, for the function y=x 
are: (-4,-4), (0,0), (3,3). Those for the function y*0 are: 
(-3.464,0), (0,0), (3.464,0) while those for the function 
y=12-x2 are: (-4,-4), (-3.464,0), (3,3), (3.464,0). Each of 
the functions has more than two endpoints and they intersect 
to form three separate regions (see Figure 1 on page 5). 

If however the user indicates that the region is 
restricted to the right half plane the program discards those 
points of intersection which are in the left half plane. 
Similarly if the user indicates that the region is restricted 
to the left half plane the program discards those points of 
intersection which are in the right half plane. 
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Once the subprogram knows all the curves and their 
endpoints, Nathematiea's Plot function samples the functions 
and determines the points which define each curve. Each curve 
and each line segment is then represented by a set of {x,y} 
points. While only the endpoints are needed here, the list of 
{x,y} points are formed to be used later in subpi^gram 
inverse. The x coordinates of the points of intersection for 
all the curves are taken to form a sorted list of points. 
Each adjacent pair of points in this list are the limits of 
integration for a subregion of the total region formed by the 
intersection of the functions. 

Subprogram cart can find the area of a region bounded by 
two curves that intersect at two different points. In this 
case the user enters the two functions of the independent 
variable x and then indicates that no limits are known. An 
example of this type of problem is the region bounded by the 
two curves y=x^ + 1 and y=x^ + x (see Figure 3 on page 27) . 
Subprogram cart can find the area of the region bounded by two 
curves that intersect at one point and by a line segment 
perpendicular to the x-axis. In this case the user enters the 
two functions of the independent variable x, indicates that a 
limit is known and then enters that limit. An example of this 
type of problem is the region bounded by the two curves y=x^ 
and y=0, which intersect at the origin, and the line x=4 (see 
Figure 4 on page 27). The third case is when both the right 
and left limits are known. Here the user enters the two 
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functions of the Independent variable x, indicates that the 
limits are known and then inputs both limits. An example of 
this type of problem is when the two curves are y=l and y=l-x'^ 
with the given limits x-1 and x=4 (see Figure 5 on page 28). 


y 




Figure 4 
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Figur* 5 


The maximum number of functions of x that the user can 
input is three. Then the maximum number of subregions which 
can be formed is two. Hence the program is able to solve for 
the area of a region which has changing boundaries. For 
example the functions y*5x, y=*x and y®cos x intersect to form 
a region which consists of two subregions (see Figure 6 
below). The first subregion has y=5x as its upper curve and 
y=x as its lower curve. The second subregion has y=cos x as 
its upper curve and y=x as its lower curve. 
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Figure 6 
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For each set of limit points and corresponding subregion, 
the subprogram determines the upper and lover curves. It does 
this by selecting those curves (now represented as lists of 
{x,y) points) whose largest and smallest x coordinate values 
bound the average value of that subregion's two limit points. 
The subprogram then uses Mathematioa's MIntegrate function to 
integrate the upper curve minus the lower curve over the 
limits of that svibregion. The region is then plotted. The 
output for each subregion includes: the limits of integra¬ 
tion, the upper and lower curves, and the area. The area of 
the total region is also provided. 

In order to solve a problem, the subprogram cart must be 
able to identify that region. It was mentioned previously 
that the subprogram has limitations if the functions 
intersect to form more than one region. The region of 
interest may be such that the program can not directly solve 
the problem. An example of this type of problem is again when 
the functions are y=x, y*0 and y=12-x^ (see Figure 1, page 5) . 
Suppose that the region of interest is that region which lies 
partly in the left half plane and partly in the right half 
plane. For this type of problem the user would have to plot 
the functions using the plot-only option (subprogram 
cartplot) . When the user indicates that his problem is given 
in cartesian coordinates he is given the option to just plot 
the functions. If this option is selected the functions are 
plotted and the points where the curves intersect are provided 
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as output. From the plot of the functions and the points of 
intersection the user could determine the limits for the 
region of interest. The user could then use this information 
and the program to solve the problem. 

In the example problem above the plot (see Figure 1, page 
5) shows that the upper curve for the siibregion in the left 
half plane is y=12-x^ while the lower curve is y=Q. For the 
subregion in the right half plane the upper curve is ysl2-x^ 
and the lower curve is y=x. The user could treat each 
subregion separately and use the program to solve for the 
areas individually and then add the two areas. 

The subprogram has no difficulty in solving problems in 
which the functions bound a single region which lies, either 
in part or whole, in the left half plane. Thus while 
subprogram oart has limitations the plot-only option can be 
used in conjunction with it to solve certain of these types of 
problems. 
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VZII. SUBPROORAX INVERSE 


Subprogram Invorso integrates those problems given in 
cartesian coordinates with respect to the variable y. First, 
however, the subprogram determines if it is appropriate to 
integrate a problem with respect to y. It attempts to find 
the inverse of each function of x that was input. If not all 
the inverse functions exist stibprogram inverse stops execu¬ 
tion. If on the other hand they all exist the user is asked 
if he wants to integrate the S 2 une problem with respect to y. 
Subprogram inverse uses Nathematiea's InversePunction to find 
the inverse functions. 

Subprogram inverse works similar to subprogram cart to 
solve the problem. The y-values of the points of intersection 
for the curves (found in subprogram cart), are limits of 
integration. There can exist limits, however, that are not 
any of the y-values of the points of intersection. An example 
of this type of problem is when the functions are y-12-x^ and 
y=-3+x+x^ (see Figure 7 on page 32). These functions 
intersect at the points (-3,3) and (2.5,5.75). The first 
function is the upper curve and takes on a maximum y-value of 
12. The second function is the lower curve and takes on a 
minimum y-value of -3.25. Hence there are three subregions. 
The limits of integration for the first region are y=-3.25 and 
y=3. The limits of integration for the second region are y=3 
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and y*5.75. The limits of integration for the third region 
are y“5.75 and y*12. The y-values -3.25 and 12 are limits 
while they are not any of the y-values for the points of 
intersection. 

y 



Figure 7 


The subprogram finds the minimum and maximum values of the 
functions and includes them as limit points. In subprogram 
cart all the curves are represented by lists of points. From 
these lists subprogram inverse selects the minimum and maximum 
values which y takes on in the curves. Once all the limits of 
integration are known they are placed in a sorted list. Each 
adjacent pair of points are then limits of integration of a 
subregion of the total region. 

For each subregion the subprogram determines the right and 
left curves in the same manner in which subprogram cart finds 
the upper and lower curves. Some inverse functions have a 
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positive branch and a negative branch. These functions are 
identified as having a square root term involving the variable 
X. Subprogram Inverse checks to see if either the right curve 
or left curve for a subregion has this type of inverse 
function. If an inverse function is of this type, the 
subprogram must determine whether it is the positive or the 
negative branch. 

If it is the negative branch, the subprogram calls 
routine negbranch, which transforms the inverse function. The 
term in the inverse function involving the square root is then 
replaced with the negative of that term. For exajaple consider 
the region bounded by the functions y=0, y=-4+x^ and x=-4 (see 
Figure 8 on page 34) . The inverse of y=-4+x^ is x=y y+4. 
When the region is integrated with respect to y the right 
curve is x-J y+4 and the left curve is x*-4. But the right 
curve is a negative branch so we have x=-y y+4 as the right 
curve. In the calculation of the integral the integrand is 
the right curve minus the left curve, or (-y y+4) - (-4) = 
-y y+4+4. 

If the right curve for a subregion is the same as the 
left curve, the subprogram calls another routine bothbranchas. 
This routine works similar to the routine negbranch. It 
insures that the region integrated over is the positive branch 
minus the negative branch (or twice the positive branch). For 
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Figur* 8 


example consider the region bounded by the functions y=x^ and 
y=4 (see Figure 9 on page 35). The inverse function of y=x^ 
is x=y^. But in this problem the right curve is the positive 
branch and the left curve is the negative branch. Therefore 
the integrand is /~x - = 2 /~lc, which is the right 

curve minus the left curve. 

The area of the subregions are then found using 
Nathematiea's Nlntagrata function. The output for each 
subregion includes: the limits of integration, the left and 
right curves, and the area. The area of the total region is 
also provided. 

Subprogram Inverse as mentioned previously works only when 
all inverse functions exist. Also stated was that it only 
integrates with respect to the variable y when it is 
appropriate to do so. There are functions whose inverse 
exists but which are complicated functions to integrate 
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Figure 9 

(Matheaatiea has no problem but the calculus student might). 
For those problems involving transcendental functions the 
subprogram will not integrate with respect to y, even if such 
integration is possible. Integration with respect to the 
variable y is supposed to save the calculus student time. 
Usually the inverse of a function containing transcendental 
functions is complicated and more difficult to integrate. 





IX. SUBPROGRAM CARTPLOT 


Subprogram eartplot is an option which allows the user to 
graph functions given in cartesian coordinates. After the 
user indicates that a problem is given in cartesian 
coordinates, subprogram graph asks if the user just wants to 
plot the functions. If the user selects this option sub¬ 
program graph calls up siibprogram eartplot. The user first 
enters all functions of the independent variable x and then 
those lines perpendicular to the x-axis. The user enters the 
ranges for both the independent variable x, and the dependent 
variable y. The functions are then plotted and the points 
where the curves intersect are given. 

This plot-only option can be used to estimate the coordin¬ 
ates of where the functions intersect and thus the limits of 
integration. There is no limit on the number of functions 
which can be plotted on one graph with this option. A single 
function can be plotted if desired. Subprogram eartplot uses 
Mathematlca's Plot function to plot all graphs. 

Subprogram eartplot can be used by itself without invoking 
the main program. This is valuable in saving time if the user 
just wants to plot several graphs. In order to use this 
stand-alone option the user types eartplot and then activates 
the cell. The tutorial covers the many ways in which the cell 
can be activated. 
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Regardless how subprogram oartplot is called the first 
prompt the user sees is: 

To just plot ourwes type j, otherwise type n. 

Suppose the user vdintB to plot the functions y»x, y^o and 
y=12-x^. The user types y to select the plot-only option. 
The next message the user sees is: 

Enter the curves« one at a time. When finished type n. 

y * 

The user then types z and presses Enter to input the first 
function. The above message appears again and the user then 
types 0 and presses Enter to input the second function. Once 
again the above message appears and the user inputs the last 
function by typing 12-z^ and pressing Enter. The above 
message appears again and this time the user types n and 
presses Enter to indicate that he is finished. 

The user then sees the message: 

Enter lines perpendicular to the z-azis. When finished 
type n. 

Since there are none the user types n. If there were lines to 
enter they would have been entered one at a time as the 
functions of the independent variable x were entered above. 
The user then sees the message: 

Enter lower limit for z-range of plot. 

The user must guess what this value is as well as the upper 
limit for the x-range and both the lower and upper limits for 
the y-range. Suppose the user knows these values. For 
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example, to enter the lower limit for the x-range he then 
types -5 and presses Enter. The user then sees the message: 

Enter upper limit for z-range of plot. 

He types 5 and presses Enter and sees the message: 

Enter lover limit for y-range of plot. 

He types >5 and presses Enter and sees the message: 

Enter upper limit for y-range of plot. 

He types 12 and the plot is displayed along with the points 
where the curves intersect (see Figure 10 below). After the 
plot is displayed the user may want to change the x or y range 
of the plot. To do this the user must call up subprogram 
cartplot again and enter all the information. 
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The points where the curves intersect are: 

{{-4.,-4.), {-3.464,0.), (0.,0.), (3.,3.), (3.464,0)) 

Figure 10 
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Z. 8UBPSOGSZ1I FOLXZl 


Subprogram polarl finds the area bounded by a polar curve 
and the origin. Subprogram graph calls up subprogram polarl 
when the user indicates the problem is to find the area 
between the origin and a polar curve. Subprogram polarl can 
be used by itself without invoking the main program. In order 
to use this stand-alone option the user types polarl and then 
activates the cell. 

Regardless how subprogram polar is called, the user first 
enters the polar curve r, which bounds the region. The polar 
curve must be a function of the independent variable Q 
representing the polar angle theta. The subprogram asks the 
user if the limits of integration are known and then to enter 
them, if they are known. For those problems in which the 
limits are known in advance the subprogram then graphs the 
region and solves the problem. The subprogram uses Mathe¬ 
matical s PolarPlot function to produce all plots. The output 
includes the limits of integration, the polar curve and the 
area bounded by the curve. 

A polar curve may be negative for certain ranges of the 
independent variable Q. For example the polar curve r=sin3Q 
takes on negative values in the range Q=ir/3 to Q=25r/3. A 
polar curve may also be undefined for certain ranges of the 
independent variable Q (i.e., when equal to the square root of 
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a negative number). The polar curve r^J cos 2Q is defined 
only for Q between Q=-ir/4 and Q=ir/4 and for Q between Q=3ir/4 
and Q=5ir/4 (see Figure 11 below). 
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Figure ll 

The subprogram determines for which values of Q the 
function r is equal to zero. It uses Mathenatica's NSolve 
function to do this. For all regions between these values of 
Q (where r=0) and the values of Q entered as limits, the 
function r must be positive, negative, or undefined. The 
subprogram determines the case for each such region. 

When the function r is defined for a region the subprogram 
uses Mathematica's NIntagrate function to find the area of 
that-region. The limits of integration are those values of Q 
which define the region. Since the problem is to find the 
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area between the polar curve and the origin the absolute value 
of the evaluated integral is taken. Therefore a region where 
the polar curve r is negative contributes positive area to the 
total area between the curve and the origin. If the function 
r is not defined over a region, no area is contributed. 

If the user does not know the limits of integration the 
polar curve is plotted for Q=0 to Q»2ir. The plot of the curve 
is displayed along with those Q values for which the function 
r is equal to zero. After a brief delay the user is told that 
the polar curve is plotted from Q^O to Qs2ir and is asked to 
input the limits of integration. From this point on the 
subprogram works exactly like it did above when the limits of 
integration were known in advance. 

If the area bounded by a polar curve and the origin is 
made up of more than one region, the limits of integration for 
all those regions are provided as output. For example the 
area between the polar curve r-J cos 2Q and the origin for Q=0 
to Q=3jr/2 consists of two separate regions (see Figure 12 on 
page 42). The limits of integration for the first region are 
al=0. and bl=ir/4 while the limits of integration for the 
second region are a2=3ir/4 and b2=5r/4. These limits are 
displayed along with the graph of the region and the total 
area. 
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•The subprogram does not tell how much area is contributed 
by each region. If the user wants to know this, he could 
treat each region as a separate problem. The program could 
then be used to determine the areas individually, one at a 
time. 
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ZZ. SUBPROGRAM POLPLOTl 


Subprogrui polplotl is an option which allows the user to 
graph the area between the origin and a polar curve. After 
the user indicates that a problem is given in polar 
coordinates, subprograa graph asks if the user just wants to 
plot the functions. If the user selects this plot-only option 
the program then asks the user to select the type of plot. 
The choices are: (1) the plot of the area between the origin 
and a polar curve, and (2) the plot of the area between two 
polar curves. When the user selects (1), subprogram graph 
calls up svibprograB polplotl. 

The user enters the polar curve which bounds the region 
and then indicates if the limits are known. If the limits are 
not known, the polar curve is plotted for the polar angle Q=0 
to Q=2n. Those points where the polar curve is equal to zero 
are displayed along with the graph. The user is then given an 
option to change the limits the user selects this option 
the polar curve is then replotted for the new limits of the 
polar angle Q, which the user enters. Subprogram polplotl 
uses Mathematica's PolarPlot function to plot the graphs. 

Subprogram polplotl can be used by itself without invoking 
the main program. In order to use this stand-alone option the 
user types polplotl and then activates the cell. 
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Regardless how aiobprograa polplotl is called the first 

prompt the user sees is: 

Bntsr the polar curve which bounds the region. 

If the user wanted to plot the polar curve r= 1 ~ cos Q, he 
would type 1 - eos[Q] and press Enter. The next message the 
user would see is: 

If the limits are known type y, otherwise type n. 

Suppose the user types n. The polar curve is then displayed 
along with those points where the polar curve is equal to zero 
(see Figure 13 below). 


y 



The polar curve r is equal to zero when Q is; (0, 2 Pi) 


Figure 13 
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After a brief delay the user sees the message: 

The polar curve is plotted from Q=o to Q=2 Pi. If you 
want to change the limits type y, otherwise type n. 


Now suppose the user wants to change the lower limit to Q=k/2 
and the upper limit to Q=3ir/2. He types y and the next 
message that appears is: 

Input new lower limit. 

The user types pi/2 and presses Enter. The next message 
that appears is: 

Input new upper limit. 

The user types 3 pi/2 and presses Enter. The plot of the 
polar curve is then displayed for the new range of the polar 
angle Q (see Figure 14 below). 
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ZII. SUBPROGRAM POLAR2 


Subprograa polar2 finds the area between two polar curves. 
Subprograa graph calls up subprograa polar2 when the user 
indicates the problem is to find the area bounded between two 
polar curves. Subprogram polar2 can be used by itself without 
invoking the main program. In order to use this stand-alone 
option the user types polar2 and then activates the cell. 

Regardless how subprogram polar2 is called, the subprogram 
directs the user to enter the two polar curves. These polar 
curves must be functions of the independent variable Q. After 
the two polar curves have been entered the user is asked if 
the limits of integration are known and to input them if 
known. 

The program determines those values of Q where the curves 
intersect. The functions which define the curves are 
represented as Taylor Series expansions. The subprogram uses 
Mathematicals MRoots function to obtain the approximate values 
of Q for which the curves intersect. These values are then 
used as starting points in Matbematlca's FindRoot function. 
This function uses Newton's Method to find the points of 
intersection. 

If the limits are not known the subprogram then plots the 
two polar curves for Q=0 to Q=2k. All plots are produced 
using Mathematiea's PolarPlot function. The values of Q for 
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which the two polar curves intersect are also provided. 
After a brief delay the user is told that the curves are 
plotted for Q=0 to Q-2ir. The subprogram asks the user if he 
wants to change these limits. If the user desires he can then 
input the new limits. If the limits of integration are known 
and inputted initially, the program plots the two polar curves 
over this region. 

Once the limits are entered the subprogram calls up a 
number of routines which check for certain types of 
conditions. Routine negcheck determines if the polar curves 
are negative. The user is told if one or more of the curves 
are negative in the range Q=lower limit to Q=upper limit. The 
user is also told that the program does not handle this type 
of problem. Routine complezcheek determines the ranges of Q 
for which the curves are undefined (if they are). Routine 
diffcheck checks to see if the outer and inner polar curve 
changes between the upper and lower limits. If the outer 
curve and inner curve changes the subprogram lets the user 
know this. The user is told that a point of intersection 
cannot lie between the two limits. The user is then directed 
to input new limits. 

Once all conditions have been checked and the program 
determines that the problem is well defined it calculates the 
area between the two polar curves. The subprogram has two 
routines which it uses to find the area. If both curves are 
defined over the entire interval routine findarea2 is used. 
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This routine uses NatheBStica's NIntegrata function to 
evaluate the single integral. If hivever the inner curve is 
undefined over certain ranges of Q the routine findareal is 
used to find the area. 

An example is when the outer curve is the circle r=l and 
the inner curve is r=y cos 2Q (see Figure 15 below). Here the 
inner curve is not defined for the polar angle Q=ir/4 to Q=3ir/4 
and for Q=-3ir/4 to -ir/4. Routine findareal first finds the 
area between the outer polar curve and the origin and then 
finds the area between the inner polar curve and the origin. 
The difference between these two areas is the area between the 
two polar curves. 
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Figure 15 
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The region bounded by the two polar curves is then 
plotted. The output includes the limits of integration, the 
outer curve, the inner curve, and the total area of the 
region. The program does not provide the areas of any 
subregions. If the user wanted to know this, each subregion 
could be treated as a separate problem. The program could 
then be used to find the area of each individual subregion, 
one at a time. 
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XIII. SUBPROGRAM POLPLOT2 


Subprograa polplot2 is an option which allows the user to 
graph the area between two polar curves. After the user 
indicates that a problem is given in polar coordinates, 
subprogram graph asks if the user just wants to plot the 
functions. If the user selects this plot-only option the 
program then asks the user to select the type of plot. The 
choices are: (1) the plot of the area between the origin and 
a polar curve, and (2) the plot of the area between two polar 
curves. When the user selects (2), subprogram graph calls up 
subprogram polplot2. 

The user enters the two polar curves which bound the 
region and then indicates if the limits are known. If the 
limits are not known, the polar curves are plotted for the 
polar angle Q=0 to Q=2ir. Those points where the two polar 
curves intersect are displayed along with the graph. The user 
is then given an option to change the limits. If the user 
selects this option the polar curves are then replotted for 
the new limits of the polar angle Q, which the user enters. 
Subprogram polplot2 uses Mathematicals PolarPlot function to 
plot the graphs. 

Sxibprogram polplot2 can be used by itself without invoking 
the main program. In order to use this stand-alone option the 
user types polplot2 and then activates the cell. 


50 





Regardless how subprograa polplot2 is called the first 
prompt the user sees is: 

Enter one of the polar curves which bounds the region. 

Suppose the user wanted to plot the area between the two polar 
curves r= 2 and r= 1 + sin Q. To enter the first polar curve 
the user would type 2 and press Enter. The next message the 
user would see is: 

Enter the other polar curves which bounds the region. 

The user would type 1 •t- sin[Q] and press Enter to enter the 
second polar curve. The next message the user would see is: 

If the limits are known type j, otherwise type n. 

Suppose the user types n. The polar curves are then displayed 
along with those points where the polar curves intersect (see 
Figure 16 on page 52). 

After a brief delay the user sees the message: 

The polar curves are plotted from Q=o to Q=2 Pi. If you 
want to change the limits type y, otherwise type n. 

Now suppose the user wants to change the upper limit to 
Q=ir. He types y and the next message that appears is: 

Input new lower limit. 

The user types 0 and presses Enter. The next message that 
appears is: 

Input new upper limit. 

The user types pi and presses Enter. The plot of the two 
polar curves is then displayed for the new range of the 
polar angle Q (see Figure 17 on page 52). 
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The two polar curves intersect at Q 


Figure 16 

y 



Figure 17 


= Pi 
~2 
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ZIV. SUBPROGRAM EXAMPLES 


The file named examples.ma contains the Matbematica code 
for each of the graphs found in the lab book. All of the 
graphs are in a package and can be displayed by the subprogram 
examples. When in the main program the first prompt says: 

To display a graph that is in the lab book type the 
example number, otherwise type n. 

The user can then type a number from 1 to 23 to call up 
subprogram examples. The subprogram then displays the graph 
which corresponds to that number example problem. 

There may be only limited value in displaying two- 
dimensional graphs which are already in the lab book. 
However, when this integral package is expanded to include the 
three-dimensional case this could be a very useful option. 
Mathsmatica has a sophisticated three-dimensional graphing 
capability which uses a three-color lighting system to 
illuminate three dimensional surfaces. Matbematica also has 
the ability to change the point of view interactively to 
explore three-dimensional graphics. For these types of 
problems being able to display the graphs becomes important. 
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ZV. rUKCTIOllS AND NOTATION 


A. FUNCTIONS 

When inputting information both the program and 
Methematice are case sensitive. Some transcendental functions 
are defined within the program code and must be in lower case 
letters. These transcendental functions are represented as 
Taylor Series expansions. They were needed to solve algebrai¬ 
cally for where two curves intersect, using NRoots. These 
approximate points were then used as starting points for 
Matheaatioa's FlndRoot function which uses Newton's Method. 

Mathematlea has its own built-in transcendental functions 
which in most instances are spelled exactly the same. The 
only difference is that all Matheaatica functions begin with 
an upper case letter. 

Below are those functions which the program recognizes: 


sin[x] 

sec[x] 

cos[x] 

exp[x] 

tan[x] 

ln[x] 

cot[x] 

sqrt[x] 

csc[x] 



B. NOTATION 

1. Brackets and Parentheses 

Square brackets and parentheses are intended for 
different purposes. Square brackets are used for specifying 





ar 9 \iments of functions. Parentheses are used for grouping. 
Without parentheses, multiplication and division have a higher 
precedence than addition and subtraction. 

2. Kathematioal Symbols 

The standard mathematical operations are referred to 
with symbols. Those symbols used are: 


Mathematical Symbol Operation 


+ plus, add 

minus, subtract 
* times, multiply 

/ divide 

^ power 


3. Other Symbols 

The symbol Q is used for theta in problems given in 
polar coordinates. All trigonometric functions must be 
functions of the variable Q. The program does not recognize 
any other such variable. The symbol pi represents the 
numerical value of the mathematical constant n. 
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XVI. rUTURB WORK 


Right now the program just handles two-dimensional 
problems. The package could be extended to include three- 
dimensional problems. The three-dimensional case would 
include problems given in cartesian as well as cylindrical and 
spherical coordinate. Three-dimensional problems could also 
be included in the lab book. Subprogram examples then could 
display the graphs for these example problems. Mathematicals 
representation of three-dimensional plots would give the 
student a feel for the depth of the region and insight into 
how the surfaces intersect. 

Currently there is no flexibility in defining the 
dependent and independent variables for either cartesian or 
polar coordinate problems. In cartesian coordinates, 
functions are of the independent variable x and the dependent 
variable is y. In polar coordinates the polar function r must 
always be a function of the independent variable Q. The 
program could be improved by allowing the user to determine 
the variable names (i.e., the user may want the dependent 
variable to be t). 

The focus of all future work should support the intent of 
the integral package, which is to teach the calculus student 
how to find the limits of integration for integral problems. 
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As such, all future work should be to improve upon the 
package's ability to do this. 
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APPENDIX A 


LAB BOOK 


I. INTRODUCTION 

One of the most difficult topics in calculus is deciding 
on the limits of integration. This lab book shows the 
calculus student how a given domain is plotted and teaches the 
student how to find the limits of integration when evaluating 
two-dimensional Integrals. 

A. Integration 

Integral calculus is the mathematics we use to find 
lengths, areas, and volumes of irregular shapes and to 
calculate the average values of functions. The development of 
integral calculus starts from the calculation of areas. 


y 



We can use the integral calculus to find the areas of 
regions like the shaded one here. 
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B. Lab Book organitation 

This package uses numerous example problems for the 
two dimensional case. It considers cartesian as well as 
polar coordinates. Almost all of the problems begin with a 
description of a given domain over which an integral is to be 
evaluated. The given domain is plotted and the way in which 
the limits of integration are found is shown. Finally the 
integrals are set up and evaluated. 


II. CAKTB8IAM COORDIMATB SYSTEM 

A. Finding the Area under the Graph of a Monnnegative 
Continuous Function. 

How to Find the Area under the Graph of a 
Non-negative Continuous Function y=f(x) from a to 
b. 

STBP 1: Find an antiderivative F(x) of f(x). 

STEP 2: Calculate F(b)-F(a). This number will be 
the area under the curve from a to b. 

Example 1 

Find the area under the curve y-x^ from x=0 to 
x=»l. 


y 
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Solution. In this example f(x)»x^, a*0, and b=l. 

He start with a graph of the function and identify 
the area of the region that lies between the curve 
and the x-axis from a to b. We then find the area 
in two steps. 


STEP 1: Find an antiderivative F(x) of f(x)=x*. 


x^ 

F(X) - - 
3 

STEP 2: Calculate F(1)-F(0). 


( 1 )’ ( 0 )’ 1 1 

F(1)-F(0) =-—--0 = - 

3 3 3 3 


The area is 1/3. 

Example 2 

Find the area under one arch of the curve y=cos x. 


y 



Solution. In this example f(x)=cos x, however a and 
b are not explictly given. We graph the function 
and observe that it crosses the x-axis at x=-ir/2 and 
x=ir/2. Note that cos{-rr/2) =0 and cos(r/2) =0. 
Thus we have a--ir/2 and b=x/2. We then find the 
area in two steps. 
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8TBP 1: Find an antiderivative F(x) of f(x)«cos x. 
F(x)=ssin X 

BTBP 2: Calculate F(-jr/2)-F(»/2) . 

F(-ir/2) ~F(ir/2) =sin(-x/2) -sin(ir/2) =1- (-1)=2 
The area is 2. 

B. The Average value of a Function 

The average value of the function f on [a,b] is the 
integral of f divided by the length of the integral. 


Average value of f on [a,b] is - J f(x)dx. 

b-a a 


If f is continuous and nonnegative on [a,b], its 
average value is the height of a rectangle whose 
area. 


f(c)(b-a)= 


Jf(x)dx, 
a 


is the area under the graph of f from a to b. 

y 
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Bzupl* 3 


Find the average value of f(x)= ^4 - , on 

the on the interval [-2, 2]. 

Solution: We graph the integrand f(x) over the 
interval of integration [-2, 2] and see that the 
graph is a semicircle of radius 2. The area between 
the semicircle and the x-axis is 


1 1 

Area = (-)jr r^ (-)»(2)*=2ir. 

2 2 

Because the area is also the value of the Integral 
of f from x=-2 to x=2, the average value of f on 
t-2, 2] is 


Average value =■ 


2 -(- 2 ) 


I \/4 - x^ dx=(-) 2w= 

-2 4 2 


y 2 

y * 4 - X 



Example 4 

Find the average value of f(x)=2x on the interval 
[1,4]. 

Solution. We graph the integrand f(x) over the 
interval of integration [1,4] and we see that the 
graph of the area is a trapezoid whose base is on 
the x-axis. The area of a trapezoid is 
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1 1 

Area® - base(hl + h2)= -(3)(8 + 2)®15. 

2 2 

Because the area is also the value of the integral 
of f from x=l to x=4, the average value of f on [1, 
4] is 



Area Between Curves 

This section shows how to find the area of a region 
in the coordinate plane by integrating the functions 
that define the region's boundaries. 

Definition 

If functions f and g are continuous and if f(x) is 
greater than or equal to g(x) throughout the 
interval [a,b] such that x is in [a,b], then the 
area of the region between the curves y=f(x) and 
y=g(x) from a to b is the integral of (f-g) from a 
to b: 
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Area= J [f(x)-g(x)]dx. 
a 

How to find th« Aroa botwaon Two Curvaa 

1. Graph the curves together. This tells you which is f, 
(upper curve) and which is g, (lower curve). It also 
helps to find the limits of integration if you do not 
already know them. 

2. Find the limits of integration. 

3. Write a formula for f(x)-g(x). Simplify it if you 
can. 

4. Integrate f(x)-g(x) from a to b. The number you get is 
the area. 


Bxaapla 5 

Find the area between the curve y=cos x and the 
curve y=-sin x from x=0 to x=x/2. 

Solution. 

STEP 1; The graphs. We graph the curves together. 
The upper curve is y*cos x, so we take f(x)=cos x 
in the area formula. The lower curve is y=-sin x, 
so g(x)=“Sin x. 


y 
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8TBP 2: The limits of integration. They are 
given: 

a^O and b*ir/2. 

STBP 3: The formula for f(x)-g(x). From Step l, 
f(X) -g(x)=cos X - (-sin x) 

*cos X + sin X 

8TBP 4: Integrate f(x)-g(x) from a=0 to b=ir/2. 

ir/2 w/2 

J (cos X + sin x)dx =[sin x - cos x] 

0 0 

=[(l-0)-(0-l)]=2 

The area between the curves is 2. 

Bxample 6 

Find the area between the curves y=l and y=l-x‘^ 
for x=l to x*4. 

Solution. 

STBP 1: The graphs. We graph the curves together. 
The upper curve is y®!, so f(x)=l. The lower curve 

is y=l-x'^, so g(x)=l-x'^. 


y 
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STEP 2: The limits of integration. The limits of 
integration are given: a-1 and bs4. 

STEP 3: The formula for f(x)-g(x). 

f(x)-g(x)=l-(l-x-2)»x-2 
STEP 4: Integrate. 

Area= f (x)-g (x) ]dx= |x‘^ dx= - - 

a lx 

The area of the region is 3/4. 

D. Curves that Cross 

When a region is determined by curves that cross, 
the crossing points give the limits of integration. 

Example 7 

Find the area of the region enclosed by the 
parabola y=2-x^ and the line y=-x. 

Solution. 

STEP 1: The graphs. We graph the curves together. 

The upper curve is y=2-x^, so f{x)=2-x^. The lower 
curve is y=“X, so g(x)=-x. The x-coordinates of the 
points where the parabola and line cross are the 
limits of integration. We find them in Step 2. 

y 
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STEP 2: The limits of integration. We find the 

limits of integration by solving the equations 
y*'2-x^ and y=-x simultaneously for x: 


2-x^*-x (Equate f(x) and g(x)) 

x^-x-2=0 (Transpose) 

(x+l)(x-2)=0 (Factor) 

x=-l, x=2. (Solve) 


The region runs from x*-! on the left to x«2 on the 
right. The limits of integration are a»-l and b^2. 

STEP 3: The formula for f(x)-g(x). 

f (x)-g(x) = (2-x^)-(-x)=2+x-x2 

STEP 4: Integrate. 


Area= J [f (x)-g(x) ]dx = J (2+x-x^)dx 
a -1 


x^ x^ 2 
2x + - - - 

2 3-1 


4 8 11 

= (4 +-) - (-2 + - + -) 

2 3 2 3 


3 9 9 

=6 +-- - 

2 3 2 

The area of the region is 9/2. 
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BXUlpl* • 

Find the area of the reqion enclosed by the curve 
and the curve y=x^+x. 

Solution. 

STEP 1: The graphs. We graph the curves together. 

The upper curve is y=x*+l, so f(x)*x^+l. The lower 
curve is y=x^+x, so g(x)=x^+x. The x-coordinates 
of the points where the curves cross are the limits 
of integration. We find them in Step 2. 


y 



STEP 2; The limits of integration. We find the 
limits by solving the equations y=x’+l and y=x +x 


simultaneously for x: 


x’+l =x^+x 

(Equate f(x) and g(x)) 

x^-x^-x+1 =0 

(Transpose) 

(X+1) (x-l)2 =0 

(Factor) 

X 

II 

X 

II 

• 

(Solve) 


The region runs from x=-l on the left to x=:l on the 
right. The limits of integration are a=-l and b=l. 
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8TBP 3: The formula for f(x)-g(x). 

f (X) -g (X) =x^+l“(x^+x) *x^-x^-x+l 
STEP 4: Integrate. 

Area= f (x)-g(x) ]dx= |[x^-x^-x+l]dx 
a “1 


111 
-x‘--x^--x^+l 
4 3 2 


1 

-1 


111 , 

-(!)" - -(1)3 - -(1)2+ ij 

4 3 2 


111 
-(- 1 )‘ - -(- 1)3 - -(- 1)2 +1 

4 3 2 


111 

- +1 

4 3 2 


111 

- +-+ 1 

4 3 2 


4 

3 


The area of the region is 4/3 


Example 9 

Find the area of the region enclosed by the 
parabola y= 12 -x 2 parabola y=4+x2 . 

Solution. 

STEP 1: The graphs. We graph the curves together. 
The upper curve is y=12-x2, so f(x)=12-x2. The lower 
curve is y=4+x2, so g(x)=4+x2. The x-coordin-ates of 
the points where the two parabolas cross are the 
limits of integration. We find them in Step 2. 


69 








8TKP 2: The limits of integration. We find the 

limits by solving the equations y=12-x^ and y=4+x^ 
simultaneously for x; 

4+x^ = 12-x^ (Equate f(x) and g(x)) 

2x^ = 8 (Transpose) 

x^ = 4 (Divide by 2) 

x=2, x=-2. (Solve) 

The region runs from x=-2 on the left to x=2 on the 
right. The limits of integration are a=-2 and b=2. 

STEP 3: The formula for f(x)-g(x). 

f (x)-g(x)=12-x^-(4+x^)=8-2x^ 

STEP 4: Integrate. 

Area= 5^[f (x)-g(x) )dx = | [8'-2x^]dx 
a -2 
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2 2 

= 8x - - x’ 

3 -2 

2 2 62 

=[ 8 ( 2 ) - -( 2 )’] - [ 8 (- 2 ) - -(- 2 )’]= — 

3 3 3 


The are’^ of the region is 62/3. 

B. Boundaries with Changing Fomulas 

If the formula for one of the bounding curves 
changes at some point across the region, you may 
have to add two or more Integrals to find the area. 

Example 10 

Find the area of the region in the first quadrant 
bounded above by the curve y=/~x and bounded below 
by the x-axis and the line y=x-2. 

Solution. 

STEP 1: The graphs. We graph the curves together. 
The entire upp er boundary of the region consists of 
the curve y=/^ x, so f(x)=y x. The lower boundary 
consists of two curves, first y=0 for x=0 to x=2 and 
then y=x-2 for x=2 to x=4. Hence the formula for 
g(x) changes from g(x)=0 for x=0 to x=2, to g(x)=x--2 
for x=2 to x=4. 
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STEP 2: The limits of integration. The limits of 
integration for the pair f{x)=/~x and g(x)=0 are a=0 
and b=2. For the pair t{x)=y~lic and g(x)*x-2, the 
left-hand limit is a=2 and the right- hand limit is 
the x-coordinate of the upper point where the line 
crosses the parabola. To find it, we solve the 
equations and y=x-2 simultaneously for x: 


y^=x-2 (Equate f(x) and g(x)) 

x=(x-2)^ =x^-4x+4 (Square) 

x^-5x+4=0 (Transpose) 

(x-1)(x-4)=0 (Factor) 

x=l, x=4. (Solve) 

The value x=l does not satisfy the equation y~x=x-2. 
It is an extraneous root introduced by squaring. 
The x=4 gives our upper limit of integration. 

STEP 3; The formulas for f(x)-g(x). 

For x=0 to x=2: f (x)-g(x)=y^-0=y^, 

For x=2 to x*4: f (x)-g(x)=/^-(x-2) 

=y~x-x+2. 

STEP 4; Integrate. We have two integrals to 
evaluate. Their sum is the area. 


Area= J[f(x)-g(x)]dx 
a 


J /~x dx + J[/^-x+2 3dx 
0 2 


2 3/2 2 2 3/2 X^ 4 

- - X +-X --+2X 

3 0 3 2 2 
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2 3/2 

2 3/2 (4)^ 1 


12 3/2 (2)2 1 

» -(2) + 

-(4) 

- - + 2(4) 

- 

-(2) 

- - + 2(2) 

3 

3 

2 


3 

2 


10 

3 

The area of the region is 10/3. 

Bxaaple 11 

Find the area of the region in the first quadrant 
bounded above by the curves y=x and y=2-x^ and below 
by the x-axis. 

Solution. 

STEP 1: The graphs. We graph the curves together. 
The upper boundary consists of two curves, first 

y=*x for x=0 to x=l and then y=2-x^ for x=l to 
x=y 5. Hence the formula for f(x) changes from 
f(x)=xfor x=0 to x=l, to f(x)=2-x^ for x=l to 
The entire lower boundary of the region consists of 

the curve y®0, so g(x)=o. 
y 



STEP 2: The limits of integration. For the pair 
f{x)=x and g(x)=0 the left-hand limit is a=0. The 
right-hand limit is the x-coordinate of the upper 
point where the line crosses the parabola. To find 
it, we solve the equations y=x and y=2-x^ simul¬ 
taneously for x: 
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X=2-X^ 


(Equate f(x) and 9 (x)) 

x^+x-2=0 (Transpose) 

(x+2)(x-l)=0 (Factor) 

x=-2, x=l. (Solve) 

Since X is greater than or equal to zero in the 
first quadrant we must pick x^i. The limits of 
integration are a=0 and b=l. For the pair 
f(x)=2-x^ and g(x)3E0 the left-hand limit is 1 and 
the right-hand limit is obtained by solving for 
where the curve y=2-x^ crosses the x-axis: 


2-x^*0 (Equate f(x) and g(x)) 

2=x^ (Transpose) 

x=y~2. (Solve) 

This region runs from 1 on the left to J~2 on the 
right. The limits of integration are a=l and b=y~2. 


STEP 3: The formulas for f(x)-g(x). 

For x=0 to x=l: f(x)-g(x)=x-o=x, 

For x=l to x=y^; f (x)-g(x) =2-x^-0=2-x^. 

STEP 4: Integrate. We have two integrals to 
evaluate. Their sum is the area. 




Area= J[f(x)-g(x)]dx 


= } xdx + J 


[2-x2]dx 


x2 

1 

X* 

= - 

+ 

2x- 

2 

0 

3 


( 1 )^ 

- + [ 2 /^ 
2 


1 


•] - [ 2 ( 1 ) - 


( 1 )" 

- ] 
3 


= .719 
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BZUpl* 12 


Find the area bounded by the parabola -x and the 
line y»6-x. 

Solution. 

STEP 1: The graphs. We graph the curves together. 
The upper boundary consists of two curves, first 
for x=0 to x=4 and then y=6-x for x=4 to 
x=»9. Hence the formula for f(x) changes from 
f(x)*y^ for x=0 to x*4 to f(x)“6-x for x=4 to 
x=9. The entire lower boundary of the region 
consists of so g(x)=-yx. 


y 



STEP 2: The limits of integration. For the pair 
f(x)*y^ and g(x)=-/^ the left-hand limit is 0. 
The right-hand limit is the x-coordinant of 
the upper point where the line crosses the 
parabola. To find it, we solve the equations 
y=J~x and y=6-x simultaneously for x: 

7~x=6-x (Equate f(x) and g(x)) 

x=x^-12x+36 (Square) 

x^-13x+36=0 (Transpose) 

(x-4)(x-9)=0 (Factor) 

x=4, x=9. (Solve) 
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The upper point where the parabola and line cross 
is at x=4. The limits of integration are a=0 and 
b=4. For the pair f(x)=6-x and g(x)=-y~x the 
left-hand limit is x=4 and the right-hand limit is 
the x-coordinate of the lower point where the 
parabola and line cross. We found this point when 
we solved the two equations above simultaneously 
for X. The right-hand limit is x=9. The limits of 
integration are a=4 and b=9. 

STEP 3: The formulas for f(x)-g(x). 

For x=0 to x=4; f(x)-g(x)=y~x-(-y~x) 

=2y~x, 

For x=4 to x=9: f(x)-g(x)=6-x-(-7 x) 

=6-x+/~x. 

STEP 4: Integrate. We have two integrals to 
evaluate. Their sum is the area. 

Area=5^[f (X) -g(x) ]dx= ^2/Hc dx + J[6-x+y^]dx 
a 0 4 

4 3/2 4 x^ 2 3/2 9 

*-x + 6X--+-X 

3 0 2 3 4 

4 3/2 (9)2 2 3/2 (4)2 2 3/2 

=-(4) +[6(9) - - + -(9) ] - [6(4) - - + -(9) ] 
3 2 3 2 3 

=20.833 

The area of the region is 20.833. 

F. Integrating with Respect to y 

When a region's bounding curves are described by 
giving x as a function of y, the basic formula 
changes, 

For regions like these. 
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y 





use the formula: 


Area= 


J[f(y)-g(y)l^y. 

a 


The only difference is that we are now integrating with 
respect to y instead of x. We can sometimes save time by 
doing so. The basic steps are the same as before. 
Example 13 

Find the area of the region between the curves x=y^ and 
x=y+2 in the first quadrant. 

Solution. 

STEP 1: The graphs. We graph the curves together. The 
right-hand curve is x=y+2, so f(y)=y+2. The left-hand 
curve is x=y^, so g(y)=y^. 
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STEP 2: The limits of integration. The lower limit of 
integration is y=0. The upper limit is the y-coordinate 
of the upper point where the line crosses the parabola. 
We find it by solving the equations x=y+2 and x=y^ 
simultaneously for y: 

y+2=y^ (Equate f(y) and g(y)) 
y^+y+2=0 (Transpose) 


(y+1)(y-2)=0 (Factor) 


y~-l, y=2. (Solve) 


The upper limit of integration is 2. (The value y=-l 
gives the point of intersection below the x-axis.) 


STEP 3; The formula for f(y)-g(y). 

f(y)-g(y)=y+2“y^ 

STEP 4; Integrate. 

Area= J[f(y)-g(y)]dy= f[2+y-y2]dy 
a 0 


yZ y3 

2y + - - - 
2 3 


2 4 8 

= 4 +- 

0 2 3 


10 

3 


The area of the region is 10/3. 

Example 14 

Find t he ar ea of the region between the curves x=y 
and x=y i-y in the first quadrant. 

Solution. 


STEP 1: The graphs. We graph the curves together. The 
right-hand curve is x=y 2-y, so f(y)=y 2-y. The left-hand 
curve is x=y, so g(y)=y. 
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y 



X 

STEP 2: The limits of integration. The lower limit of 
integration is y=0. The upper limit is the y-coordinate 
of the upper point where the line crosses the parabola. 
We find it by solving the equations x=y and x=y i-y 
simultaneously for y: 

y=7 2-y (Equate f(y)-g(y)) 

y2=2-y (Square) 

y^+y-2=0 (Transpose) 

(y+2)(y-l)=0 (Factor) 

y=-2, y=l. (Solve) 

The upper limit of integration is l. (The value y=-l gives 
the point of intersection below the x-axis.) The limits 
of integration are a=0 and b=l. 

STEP 3: The formula for f(y)-g(y). 

f(y)-g(y)=/^^ - y 

STEP 4: Integrate. 


Area= J[f(y)-g(y)]dy = J[7”5-y - y)dy 
a 0 


2 3/2 y^ 1 

= - -(2-y) - - 

3 2 0 
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2 3/2 (1)* 212 3/2 

=t- -(2-1) - -)*- - + -(2) - 0=.719 

3 2 3 2 3 

The area of the region is .719. 

Example 15 

Find the area bounded by the parabola y*=x and the line 
y=6-x. 

Solution. 

STEP 1: The graphs. We graph the curves together. The 
right-hand curve is x=6-y, so f(y)=6-y. The left-hand 
curve is x=y^, so g(y)=y*. 


y 



STEP 2: The limits of integration. The lower limit of 
integration is the y-coordinate of the lower point where 
the parabola crosses the line. The upper limit is the 
y-coordinate of the upper point where the parabola and the 
line cross. We find these points by solving the equations 
x=y^ and x=6-y simultaneously for y: 

y^=6-y (Equate f(y) and g(y)) 

y^+y-6=0 (Transpose) 

(y+3)(y-2)=0 (Factor) 

y=-3, y=2. (Solve) 

The upper limit of integration is 2 while the lower limit 
is -3. The limits of integration are a=-3 and b=2. 

STEP 3: The formula for f(y)-g(y). 

f(y)-g(y)=6-y-y^ 
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STEP 4: Integrate. 


Area= 


J[f(y)-g(y)]dy 

a 



[6-y-y2]dy 


6y 


y2 y3 

2 3 


=[ 6 ( 2 ) 


( 2)2 ( 2)5 (- 3)2 (- 3)5 

- - - ] - [6(-3) --— ] 

2 3 2 3 


8 9 

=[12 - 2 -] - [-18 - + 9) = 20.833 

3 2 

The area bounded by the curves is 20.833. 
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III. POLAR COORDINATE SYSTEM 


This section shows how to calculate areas of plane 
regions and lengths of curves. The general methods for 
setting up the integrals are the same as for cartesian 
coordinates, although the resulting formulas are somewhat 
different. 

A. Area in the Plane 



We calculate areas of plane regions like the shaded one 
above using polar coordinates. 

Area Between the Origin and r=f(Q), Q between a and b: 


Area = - r^ dQ. 

a 2 

Example 16 

Find the area of the region enclosed by the cardioid 
r=2(1+cosQ). 

Solution. 

STEP 1; The graph. We graph the cardioid and determine 
that the radius OP sweeps out the region exactly once as 
Q runs from 0 to 2w. 
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STEP 2: The limits of integration. The limits of 

integration are a=0 and b=2». 

STEP 3: Determine the integrand. 

1 1 

- = -(4) (l-t-cosQ)^ = 2 (l+2cosQ+cos^Q) 

2 2 

1+COS2Q 

= 2+4cosQ+2[-] = 3 + 4COSQ+COS2Q 

2 

STEP 4; Integrate. 


J n sin2Q 2jr 

[3+4cosQ+cos2Q]dQ = 3Q+4sinQ+[-] 

0 2 0 

= 6ir+0+0=6ir 

The area of the region is 6ir. 

Example 17 

Find the area of the region bounded by the rays Q=-jr/4 and 
Q=jr/4, and the graph of the cardioid r=l+sinQ. 

Solution. 

STEP 1: The graph. The region is the portion of the 
cardioid swept out by the radius of length r=l+sinQ as Q 
increases from Q=-ir/4 to Q=:jr/4. 
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y 



STEP 2: The limits of integration. The limits are given; 

a=-ir/4 and b=jr/4. 

STEP 3: Determine the integrand. 

11 1 
- r^ = -[1+sinQ]^ = -[l+2sinQ+sin^Q] 

2 2 2 

1 11 

= -[l+2sinQ+(- - -COS2Q) ] 

2 2 2 

13 1 

= -[- + 2sinQ - -COS2Q] 

2 2 2 

3 1 

= - + sinQ - -COS2Q 

4 4 

STEP 4: Integrate. 


1 ,r/4 3 1 

- r^ dQ = J [- + sinQ - -cos2Q]dQ 
a 2 -ir/4 4 4 

3 1 
= -Q - cosQ - -sin2Q 

4 8 

The area of the region is .928. 

Example 18 

Find the area inside the smaller loop of the limacon 
r=2cosQ+l. 


7r/4 

= .928 

-jr/4 
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Solution. 


STEP 1; The graph. After sketching the curve we see that 
the smaller loop is traced out by the point (r,Q) as Q 
increases from Q=2ir/3 to Q=4ir/3. 

y 



STEP 2: The limits of integration. Since the curve is 
symmetric about the x-axis, we may calculate the area of 
the shaded half of the inner loop by integrating from 
Q=2ir/3 to Q=ir. The limits are then a=2)r/3 and b=ir. The 
area we seek will be twice the value of the resulting 
integral; 


Area= 2 


1 

2X/3 


- r^ dQ = / r* dQ. 
2 2»r/3 


STEP 3: Determine the integrand. 


r^ = (2cosQ +1)^ = 4cos^Q + 1 

1+COS2Q 

= 4 [ - — - ] + 4cosQ + 1 

2 

= 2 + 2cos2Q + 4cosQ + 1 
= 3 + 2cos2Q + 4cosQ 


STEP 4: Integrate. 

w 

Area= J [3+2cos2Q+4cosQ3dQ 
2x/3 
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w 


3Q+sin2Q-»-48inQ 


2ir/3 


y~3 4/-5 

= 3ir - [2ir - - + -] 

2 2 


3/~T 


= ir - 


The area of the smaller loop is ir - 3y '3/2. 

Example 19 

Find the area of the region enclosed by the graph of the 
polar equation r=sin3Q. 

Solution. 


STEP 1: The graph. When we graph the region we see that 
it is a three-leaved rose. 
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STEP 2: Limits of integration. Since the three leaves 
are congruent, and since the leaf determined by Q from Q=0 
to Q=*ir/3 is symmetric about the ray Q=n/6, we may obtain 
the area of the figure by calculating the area of the half 
leaf determined by Q=0 to Q=ir/6 and multiplying the result 
by 6. The limits of integration are a=0 and b=)r/6. 


ir/6 1 

Area= 6 / - r^ dQ 

0 2 

STEP 3: Determine the integrand. 

1 

6 [- r^] = 3r2 = 3(sin3Q)2 
2 

11 3 3 

= 3[- - -COS6Q] = - - -COS6Q 

2 2 2 2 

STEP 4: Integrate. 


ir/6 n 
0 4 

The area of the region is w/A. 

Example 20 


n/S 3 3 

Area= J - -cos6Q]dQ = 

0 2 2 


3 1 

-Q - -sin6Q 
2 4 


Find the area of the region enclosed by the lemniscate r^ 
=cos2Q. 


Solution. 


STEP 1: The graph. By extracting square roots on both 
sides of the equation we obtain r=y cos2Q and we graph the 
figure. ^ 



87 



8TBP 2: The limits of integration. We note that the 

function r is defined only for Q between Q=-w/4 and Q=ir/4 
and for Q between Q=3ir/4 and Q=-ir/4. As Q ranges through 
these two intervals the two lobes of the lemniscate are 
traced out. The rays Q=-jr/4 and Q-w/i determine half the 
area of the lemniscate. We may therefore integrate from 
Q=-ir/4 to Q=ir/4 to obtain the area of one lobe and double 
the result. The limits are then a=-ir/4 and b=ir/4. 

STEP 3: Determine the integrand. 

r^ = C082Q 

STEP 4: Integerate. 


Area= 



1 

cos2QdQ = -sin2Q 


w/4 


= 1 


-r/4 

The area of 

B. Area of 


2 


-ir/4 


the region is 1. 

Region which lies Between two Polar Curves 


y 



To find the area of a region which lies between two polar 
curves from Q=a to Q=b, we subtract the integral of 
(1/2)rl^ dQ from the integral of (1/2)r2^ dQ. This leads 
to the following formula. 

Area of the Region: rl(Q) less than or equal to r and r 
less than or equal to r2(Q), Q from Q=a to Q=b 
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b 1 b 1 b 1 

Area= J - r2^ dQ - J - rl^ dQ = J -(r2^ -rl^)dQ. 
a 2 a 2 a 2 

Bxaapl* 21 

Find the area of the region that lies inside the circle 
r=l and outside the cardiod r=l-cosQ. 

Solution. 

STEP 1: The graphs. We sketch the region to determine 
its boundaries and find the limits of integration. The 
outer curve is r2=l, and the inner curve is rl“l-cosQ. 

y 



STEP 2: Limits of integration. We see that Q runs from 
-ir/2 to ir/2. Because of symmetry we find the area for Q=0 
to Q=jr/2. The limits of integration are a=-ir/2 and b=»r/2. 


jr/2 1 ir/2 1 

Area= J -(r2^ -rl2)dQ = 2 J -(r2^ -rl^)dQ 
-ir/2 2 0 2 


k/2 

= J (r2^ - rl^)dQ 
0 

STEP 3: Determine the integrand. 


r2^-rl^ = 1-(l-2cosQ+cos^Q) = 2cosQ-cos^Q 
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COS2Q+1 


= 2cosQ - 


2 


STEP 4: Integrate. 


«r/2 COS2Q+1 

Area= J[2cosQ - -]dQ 

0 2 


2sinQ - 


sln2Q 


n/2 

0 


w 

* 2 - - 

4 


The area of the region is 2 - ir/4. 

Example 22 

Find the area of the region lying inside the circle 
r=3cosQ and outside the cardioid r=l+cosQ. 

Solution. 

STEP 1: The graphs. We sketch the graphs of the two 
equations r2=3cosQ and rl=l+cosQ. 

y 



STEP 2: The limits of integration. By solving the two 
equations simultaneously, we find the point of inter¬ 
section; 
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3cosQ«l-»-cosQ 


(Equate rl and r2) 
(Transpose) 


2cosQsl 
1 

cosQ=- (Divide) 

2 

Q=’f/3, Q=-ir/3. (Solve) 

The limits of integration are a=-)r/3 and b=ir/3. 
STEP 3: Determine the integrand. 


-(r2^ - rl2) = -((3 cosQ) 2 - (1 + cosQ)^] 

2 2 

1 

= -[9cos^Q -(1 + 2cosQ + cos^Q) ] 
2 

1 

= -[Scos^Q - 1 - 2cosQ] 

2 

1 

= 4COS^Q - - - COSQ 
2 

11 1 

= 4(- + -COS2Q] - - - cosQ 

2 2 2 

1 

= 2 + 2cos2Q - - - cosQ 
2 

3 

= - + 2COS2Q - cosQ 
2 

STEP 4: Integrate. 

7r/3 3 

Area= J [- + 2cos2Q - cosQ]dQ 
-jr/3 2 

3 ir/3 

= “Q + sin2Q - sinQ 

2 -rr/3 


= Pi 




The area of the region is w. 

C. The Length of a Curve 

We calculate the length of a curve r=f(Q), Q from a to b 
by expressing the differential ds=y~(dx^ +dy^) in terras of 
Q and integrating from a to b. 

If r=f(Q) has a continuous first derivative for Q from a 
to b and if the point P(r,Q) traces the curve r=f(Q) 
exactly once as Q runs from a to b, then the length of the 
curve is given by: 


W dr^ 

r^ + (—) dQ. 
a dQ 

Example 23 

Find the length of the cardioid r=l-cosQ. 

Solution. 

STEP 1: The graph. We sketch the cardioid to determine 
the limits of integration. The point P(r,Q) starts at the 
origin and traces the curve once, counterclockwise. 

y 

P(r.Q) 



STEP 2: Limits of integration. We see that Q runs from 
0 to 2n, so the limits of integration are a=0 and b=2x. 

STEP 3; Determine the integrand. 

dr 

r=l-cosQ, — = sinQ, 
dQ 






dr^ 

+ (—) “(1-cosQ)^ +(sinQ)^ =l-2cosQ+cos^Q+sin^Q 
dQ 

=l-2cosQ+l*2-2cosQ 
STEP 4 ; Integrate. 


Length* + (—) dQ = | 

a * dQ 0 


2-2COSQ dQ 




The length of the curve is 8. 
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APPENDIX B 


TUTORIAL 


A. INTRODUCTION 

This tutorial instructs the user on how to start Matbe- 
■atlca and how to load and use the program. The tutorial 
includes three example problems that are intended to show the 
user how the program works. The first example problem is an 
integral problem given in cartesian coordinates. In the first 
problem, the program is used to find the area of a region by 
first integrating with respect to the variable x and then with 
respect to the variable y. The second problem is an integral 
problem given in polar coordinates. In this problem the 
program finds the area between a polar curve and the origin. 
The final example problem is also given in polar coordinates. 
In this problem the program is used to find the area of the 
region bounded by two polar curves. These three example 
problems are representative problems that the program can 
solve. They show some of the program's many capabilities. 
The tutorial, while not comprehensive, gives the user some 
idea of the types of problems the program can solve. It 
demonstrates how the program prompts the user for input and 
how that input is entered. 
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B. STARTIMO MATHEMATICA AND LOADING PROGRAM 

To Start Nathaaatica, double-click the Matheaatica icon. 
Mathanatica documents are called Notebooks. When you start 
Mathematica, an empty Notebook window appears on your screen. 
As soon as you start typing, a cell is created, indicated by 
a cell bracket along the right edge of the Notebook window. 
Each piece of information in a Notebook is contained in a 
separate cell. To load the program and start the kernel first 
type: <<integral.ma. 

You now must activate the cell. There are several ways to 
do this. You can use the Mouse and select Evaluate Selection, 
in the Action menu which appears at the top of the screen (see 
Figure 18 on page 96). The three keyboard equivalents are: 

1. Shift+Enter 

2. Insert 

3. 5 on the numeric keypad. 

Once you have activated the cell you are ready to use the 
program. 
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C. BZMfPLB 1 


Find the area of the region bounded by the functions y = 
5x, y = X and y = 4 - x^ in the first quadrant. 

* Press the down arrow to create a new cell and then type 

graph. 

* For this problem press the Insert key to activate the 
cell. 

* You will see appear the message: 

To display a graph that is in the lab book type the 
example number« otherwise type n. 

Because you are going to enter our own example type n and 
press Enter. 

* You will see the message: 

Type 1 for cartesian coordinate system, or type 2 for 
polar coordinate system. 

since your first example is a problem given in cartesian 
coordinates type 1 and press Enter. 

* You will see the message: 

To just plot curves type y, otherwise type n. 

Since you do not want to just plot the functions type n. 

* You will see the message: 

Enter the functions which define the region's 
boundaries. When finished type n. 

The functions must be functions of the independent 
variable x and are entered one at a time. They can be 
inputted in any order. You are now ready to enter the 
functions. Type 5x then press Enter, type x then press 
Enter, type 4 - x^ then press Enter. 

* You will see the message: 

If any of the limits are known type y, otherwise type n. 

Since none of the limits are explicitly given you type n. 

* You will see the message: 
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If th« rsglon is rsatr^otsd to ths right half plans 
typs 7 , othsnriss typs n. 

Because we are told In the problem statement that the 
region is restricted to the first quadrant type y. 

* The region is then plotted and the program solves the 
problem. You will see the following output: 


y 



The limits of integration for the first region are x = 0 
and X = 0.702 

The upper curve for the first region is y = 5 x 

The lower curve for the first region is y = x 

The area of the first region is 0.985608 

The limits of integration for the second region are x = 

0.702 and X = 1.562 

The upper curve for the second region is y = 4 - x^ 

The lower curve for the second region is y = x 

The area of the second region is 1.31145 
The area of the total region is 2.29706 

All of the functions are invertible and and it is 
appropriate to integrate the problerit with respect to the 
variable y. As a result, after a breif delay you will 
see the message: 

Type y if you wish to integrate the same problem with 
respect to y, otherwise type n. 
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Type y and you will see the following output: 

The limits of integration for the first region are y = 0 
and y = 1.56 

The right curve for the first region is x = y 

y 

The left curve for the first region is x = - 

5 

The area of the first region is 0.973443 

The limits of integration for the second region are y = 
1.56 and y = 3.51 

The right curve for the second region is x=Sqrt[4-y] 

y 

The left curve for the second region is x = 

5 

The area of the second region is 1.32362 
The area of the total region is 2.29706 


D. EXAMPLE 2 

Find the area of one leaf o f the three-leaf rose formed by 
the polar curve r =y sin3Q. 

* Press the down arrow to create a new cell and then type 

graph. 

* Activate the cell by using the Mouse to select Evaluate 
Selection from the Action menu. 

* You will see appear the message: 

To display a graph that is in the lab booK type the 
example niimber, otherwise type n. 

Because you are going to enter our own example type n and 
press Enter. 

* You will see the message: 

Type 1 for cartesian coordinrte system/ or type 2 for 
polar coordinate system. 

Since your second example is a problem given in polar 
coordinants type 2 and press Enter. 

* You will see the message: 

To just plot curves type y, otherwise type n. 

Since you do not want to just plot the polar curve type n. 
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* You will see the message: 

Type 1 for area between the origin and a polar curve, or 
type 2 for area between two polar curves. 

Type 1 since your problem is to fii.d the area between the 
origin and the polar curve r -J sin 30 for one leaf. 

* You will see the message: 

Enter the polar curve which bounds the region. 

Type aqrt[sin[3Q]]. 

* You will see the message: 

If the limits are knotrn type y, otherwise type n. 

Since the limits are not known type n. 

* The polar curve is then plotted. You will see the 
following output: 

y 



The polar curve intersects the origin at Q: 

Pi 2 Pi 4 Pi 5 Pi 

{0, f , Pi, t t 2Pi) 

3 3 3 3 

After a brief delay you will see the message: 
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Th« polar eurvo is plotted from Q=o to Q=2 Pi. If you want 
to ohangs the limits type y, othsrwiss type n. 

Since you are only Interested in the area of one leaf 
type y. 

* You will see the message: 

Input new lower limit. 

Type 0. 

* You will see the message: 

Input new upper limit. 

Type pi/3. 

* The region in which you are interested is then plotted. 
You will see the following output: 



The limits of integration for the region which is bounded 
by the polar curve r = Sqrt[sin[3Q]] and the origin are: 


Pi 

a=0 and b=— 

3 

The area bounded by the polar curve and the origin is 
.333 
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B. EXAMPLE 3 

Find the area of the region bounded by the circle r » 1 
and the cardioid r - 1 > cos Q if the cuter curve is the 
cardioid and the inner curve is the circle. 

* Press the down arrow to create a new cell and then type 
graph. 

* Activate the cell by pressing Shift+Enter. 

* You will see appear the message: 

To display a graph that is in the lab book type the 
example number, otherwise type n. 

Because you are going to enter our own example type n and 
press Enter. 

* You will see the message: 

Type 1 for cartesian coordinate system, or type 2 for 
polar coordinate system. 

Since your third example is a problem given in polar 
coordinates type 2 and press Enter. 

* You will see the message: 

To just plot curves type y, otherwise type n. 

since you do not want to just plot the polar curves type 

n. 

* You will see the message: 

Type 1 for area between the origin and a polar curve, or 
type 2 for area between two polar curves. 

Type 2 since your problem is to find the area between two 
polar curves. 

* You will see the message: 

Enter one of the polar curves which bounds the region. 
Type 1 - cos[Q]. 
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* You will see the message: 

Enter the other polar curve which bounds the region. 

Type 1. 

* You will see the message: 

If the limits are known type y, otherwise type n. 

since we do not know where the two curves Intersect and 
what the limits of integration are type n. 

* The two polar curves will then be plotted. You will see 
the following output: 


y 



The two polar curves intersect at Q: 

Pi 3 Pi 

[ — , -) 

2 2 

* After a breif delay you will see the message: 

Input lower limit. 

Since the two curves intersect and bound two different 
regions you must be careful in selecting the correct 
limits. Type pi/2. 
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You will see the message: 
Input upper Unit. 

Type 3/2 pi. 


The region is then plotted. You will see the following 
output: 


y 



The outer polar curve is r = 1 - cos[Q] 

The inner polar curve is r = 1 

Pi 3 Pi 

The limits of integration are a=— and b= 

2 

The area of the region bounded by the two 
polar curves is 2.7854 
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APPBIIDIZ C - PROGRAM CODE 
I. FILE CART.MA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 

cart:= 

Module [{}, 

cltranscend; In [ x_] = restrictionR=n; restrictionIr=n; 
endpoints^{}; intersect={); 

inyfunc={}; rightliinit=n;leftliinit=n; type2D*l; 


While[type2D==l, { 

c2func=Input["Enter the functions which 

define the region's boundaries. When finished type n. 


y =•’]; 

If[c2func =!=n, ( 

functotal=ToExpression[StringJoin["y==","c2func"]]; 
AppendTo[itiyfunc, functotal]), type2D=0]; 

If[Length[myfunc]==3,type2 D=0] 

}]; 

myfunc=myfunc /. sqrt->Sqrt; 
myfunc=myfunc /. pi->Pi; 

If[Length[myfunc]==2, { 

limitsknown=Input["If any of the limits is known type 
y, otherwise type n."]; 

While[limitsknown ===y, 

leftlimit=Input["Input left limit if known, otherwise 
type n. 

X =" ]; 

leftlimit=leftlimit /. sqrt->Sqrt; 
leftlimit=leftlimit /. pi->Pi; 

If[leftlimit =!= n, AppendTo[myfunc, 
ToExpression[StringJoin["x==","leftlimit"]]]]; 
rightlimit=Input["Input right limit if known, otherwise 
type n. 
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X ="]; 

rightlimit=rightlimit /. sqrt->Sqrt; 
rightlimit=rightlimit /. pi->Pi; 

If trightliinit=!» n, AppendTo[iDyfunc, 
ToExpression[StringJoin["x=='*,"rightlimit*'] ] ] ]; 
limitsknovn = n] ; 

}]? 

If [leftlimit »= n && (rightlimit >0 j| rightliinit~=n), 
restrictionR=Input[''If the region is restricted to the 
right half plane type y, otherwise type n."] 

]; 

If[rightlimit === n && restrictionR===n && 

(leftlimit <0 j| leftlimit=*=n), 

restrictionL=Input[''If the region is restricted to the 
left half plane type y, otherwise type n."] 

]; 

last=Length[myfunc]; 

For[ii=l, ii<= last, 

(ival=ii; ++ii; 

For[jj=l, jj<=last, 

(jval=jj; ++jj; 

If[jvall=ival && 

(myfuncC[ival,l]]===y |j myfunc[[jval,1]]*==y), 

{ 

ln[x_]:=Log[x]; 

If[myfunc[[ival,1]]===y && myfunc[[jval,1]]===y, 

{ 

If[{myfunc[[ival]], myfunc[[jval]])=!= 

({myfunc[[ival]], myfunc[[jval]] /. Log[x_]->x}), 
{theroots={nSolve[(myfunc[[ival]],myfunc[[jval]])]) 

), 

{ 

ranscend; 

theroots=NSolve[{myfunc[[ival]], myfunc[[jval]]},{x,y}]; 
theroots=N[{x,y) /. theroots]; 
theioots=DeleteCases[DeleteCases[DeleteCases[ 
theroots, {x_Complex,y_)],{x_,y_Complex)], 
{x_Complex,y_Complex}]; 
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If[Length[theroots]=*0,{ 

theroots=Solve [ {inyfunc [ [ ival ] ], myfunc [ [ j val ] ]}, {x, y} ] ; 
theroots=N[(x,y) /. theroots]; 
theroots-DeleteCases[DeleteCases[DeleteCases[ 
theroots, {x_Complex,y_)],{x_,y^Complex)], 

{x_Complex, y_Coinplex) ] ) ] ; 

cltranscend; 

Ifttheroots=!*{x,y} && theroots 1={) && 

{myfunc[ [ival] ],myfunc[ [ jval] ] )«= 

({(myfunc[[ival]],myfunc[[jval]])) /. Sqrt[x_]->x), 

{ 

rmtranscend; 

theroots=Table[FindRoot[myfunc[[ival,2]]- 
myfunc[[jval,2]]==0, {x,theroots[[m,1]])], 

{m,Length[theroots]}]; 
theroots=x /. theroots; 

theroots=Table[{theroots[[m]] , (myfunc[[ival,2]] /. 
x-> theroots[[m]])),{m,Length[theroots]}]; 
temproots= 

Table[{TrueQ[(Abs[myfunc[[ival,2]]-myfunc[[jval,2]]] /. 

X -> theroots[[s,1]]) <=.5]},{s,Length[theroots])]; 
theroots=Transpose[{theroots,temproots}]; 

theroots=Transpose[ DeleteCases [theroots, (_, (False)), 1] ] [ [1] ] 
If[theroots== {}[[!]]/theroots={) ]; 

x=. 

)]; 

}]; 

}]; 

rmtranscend; 

If[myfunc[[ival,l]]===y && myfunc[[jval,l]]===x, 
theroots=N[{{myfunc[[jval,2]],(myfunc[[ival,2]]/. 
x->myfunc[[jval,2]])}}]]; 

If[myfunc[[ival,1]]===x && myfunc[[jval,1]]===y, 
theroots=N[{(myfunc [\ -val,2]],(myfunc[[jval,2]]/. 
x->myfunc[[ival,2]])}> J]; 

theroots=DeleteCases[DeleteCases[DeleteCases[ 
theroots, {x_Complex,y_}],{x_,y_Complex)], 

{x_Complex,y_Complex}] ; 

theroots=ReplaceAll[theroots, 0. -> 0]; 

If [ restrict ionR====y, 

theroots=Select[theroots,(#[[1]]>=0)&]]; 

If[restrictionL===y, 

theroots=Select[theroots,(#[[1]]<=0)&]]; 
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If[theroots =!= {x,y) && theroots =!= {),{ 
endpoints=Join[endpoints, theroots] }]; 

)l; 

If[leftlimit =!=n, endpoints=Select[endpoints,(#[[l]]>= 
leftlimit)i]]; 

If[rightlimit =!=n, endpoints=Select[endpoints,(#[[1]]<= 
rightlimit)&]]; 

)]; 

endpoints=N[Round[endpoints*1000]/1000]; 
endpoints=Union[endpoints,{}]; 
intersect»Join[intersect,endpoints]; 
intersect=Union[intersect,{)); 

ycoord=Sort[Union[Transpose[Union[intersect,{}]][[2]],{}]]; 
ycoord=N[Union[Round[ycoord*100]]/100]; 

Bm=Transpose[{endpoints[[!]], endpoints[[2]])]; 
mnew=Sort[mm[[1]]]; x=.; 

If[ival==l, { 

If[mnew[[l]] 1= mnew[[2]], 

{linel=Plot[myfunc[[ival,2]],{x,mnew[p]],mnew[[2]]), 
DisplayFunction-> Identity]; linell=linel[[1,1,1,1]]}, 

{1inel=Show[Graphics[{Line[endpoints])], 

DisplayFunction-> Identity]; linell=endpoints)] 

)]; 

If[ival==2, { 

If[mnew[[l]] 1= mnew[[2]], 

{line2=Plot[myfunc[[ival,2]],{x,mnew[[1]], mnew[[2]]}, 
DisplayFunction-> Identity]; Iine22=line2[[1,1,1,1]]}, 

{line2s=Show[Graphics[ {Line[endpoints]} ], 

DisplayFunction-> Identity]; line22=endpoints}] 

}]; 

If[ival==3, { 

If[mnew[[l]] != mnew[[2]], 

{line3*Plot[myfunc[[ival,2]],{x,mnew[[1]], mnew[[2]]}, 
DisplayFunction-> Identity]; Iine33=line3[[1,1,1,1]]}, 
{line3=Show[Graphics[{Line[endpoints])], 

DisplayFunction-> Identity]; line33=endpoints}] 

}]; 

If[ival==4, { 

If[mnew[[l]] != mnew[[2]], 

{line4=Plot[myfunc[[ival,2]],{x,mnew[[1]], mnew[[2]]}, 
DisplayFunction-> Identity]; Iine44=line4[[1,1,1,1]]}, 

{1ine4=Show[Graphics[{Line[endpoints]}], 

DisplayFunction-> Identity]; line44=endpoints}] 
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}]; 

If[ival**5, { 

If[mnew[[l]] !« innew[[2]], 

{line5=Plot[myfunc[ [ival,2] ], {x,iBnev[ [1] ], innew[ [2] ]), 
DisplayFunction-> Identity]; line55=line[[1,1,1,1]]), 

{line5=Show[Graphics[{Line[endpoints ]}], 

DisplayFunction-> Identity]; line55=endpoints)] 

}]; 

If[ival==6, { 

If [innew[ [1] ] iiinew[ [2] ], 

{line5=Plot[myfunc[ [ival,2] ], {x,Binew[ [1] ], i&new[ [2] ]}, 
DisplayFunction-> Identity]; Iine66=line6[[1,1,1,1]]), 

(line6=Show[Graphics[{Line[endpoints ]}], 

DisplayFunction-> Identity]; line66=endpoints)] 

}]; 

endpoints={) 

) 

]; 

If[ival==2, {Show[1inel,1ine2,DisplayFunction-> 

$DisplayFunction]; 
thelines={linell,line22))]; 

If[ival==3, {Show[1inel,1ine2,1ine3 , DisplayFunction-> 
$DisplayFunction]; 
thelines={linell,line22 , line33}}]; 

If[ival=*4, {Show[1inel,1ine2,1ine3,1ine4,DisplayFunction-> 
$DisplayFunction]; 

thelines»{linell,line22,line33,line44}}]; 

If[ival==5 , {Show[1inel , 1ine2 ,1 ine3,1ine4,1ineS, 

DisplayFunction-> $DispiayFunction]; 
thelines=(linell,line22 , line!3,line44 , lines5))]; 

If[ ival"6 , {Show[ 1 inel, line2, line3, line4, lines, line6, 
DisplayFunction-> $DisplayFunction]; 

thelines={linell,line22,line33,line44,lineSS,line66})]; 

intersect=Transpose[intersect]; 
intersect=Sort[intersect[[1]]]; 
intersect=ReplaceAll[intersect,0->0.]; 

xcoord=Union[intersect,{)]; 
xcoord=Union[xcoord/l.]; 
nntranscend; 

For[kk=l, kk<=Length[xcoord]-l, 

{kval=kk; ++kk; 

xx=(xcoord[[kval]]+xcoord[[kval+1]])/2; 

integrfunc=Select[thelines,(xx> #[[1,1]] && xx< #[[-1,1]])&]; 
placel=Position[thelines, integrfunc[[1]]][[1]][[1]]; 
place2=Position[thelines, integrfunc[[2]]][[1]][[1]]; 
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x=. ; rmtranscend; 

area=NIntegrate[Evaluate[myfunc[[placel,2]]-myfunc[[place2,2 

]]. 

{X, xcoordt [Jcval] ], xcoordt [Jcval+1] ]} ] ] ; 
cltranscend; ln[x_]».; 

If[area <0, 

{upper>myfunc[[place2]]; lower=myfunc[[placel]]), 
{upper=myfunc[[placel]]; lower=myfunc[[place2]])]; 

If[kval==l && Length[xcoord]~2, 

{Print["The limits of integration are x = ", 
xcoord[ [Icval] ]," and x = ", xcoord[ [kval+1] ]); 

Print["The upper curve is y = ", upper[[2]]]; 

Print["The lower curve is y = ", lower[[2]]]; 

Print["The area of the region is ", Abs[area]] }]; 

If[kval==l && Length[xcoord]==3, 

{Print["The limits of integration for the first region 
are x = ", xcoord[[kval]]," and x =", xcoord[[kval+1]]]; 
Print["The upper curve for the first region is y = ", 
upper[[2]]]; 

Print["The lower curve for the first region is y = ", 
lower[[2]]]; 

Print["The area of the first region is ", Abs[area]]; 
areal=Abs[area])]; 

If[kval“2 && Length[xcoord]“3, 

(Print["The limits of integration for the second region 
are x = ", xcoord[[kval]]," and x * ", xcoord[[kval+1]]]; 
Print["The upper curve for the second region is y = ", 
upper[[2]]]; 

Print["The lower curve for the second region is y = ", 
lower[[2]]]; 

Print["The area of the second region is ", Abs[area]]; 

Print["The area of the total region is ", areal+ 

Abs[area]]}]; 

)]; 

inverse; 

]; 


cltranscend;=Module[(),sin[x_]=.; cos[x_]=.; tan[x_]=.; 
cot[x_] = .; csc[x_] = .; sec[x_]='.; exp[x_] = .] 

transcend:=Module[{}, 

sin[X_]:= X - X^3/6 + X^5/120 - X^7/5040 + X^9/362880 - 
x^ll/39916800; 

COS[X_]:= 1 - X^2/2 + X^4/24 - x^6/720 + X^8/40320 - 
X^lO/3628800 + X^12/479001600 - 
X^14/87178291200 + X^16/20922789888000; 
tan[X_]:= X + X''3/3 + 2x'^5/15 + 17X''7/315 + 62x''9/2835 + 
1382x^11/155925; 

cot[x_]:= 1/X -X/3 -x^3/45 - 2x^5/945 - X^7/4725 - 
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2x^9/93555; 

csc[x_]:= 1/x + x/6 + 7x^3/360 + 31x^5/15120 + 

127x^7/604800 + 73x^9/3421440; 
sec[X_]:= 1 + X*2/2 + 5X^4/24 + 61x^6/720 + 277x^8/8064 + 
50521x^10/3628800; 

exp[X_]:= 1 + X + X^2/2 + X-3/6 + X^4/24 + X-5/120 + x-6/720 
+ 

X^7/5040 + X''8/40320 + X^9/362880] 

nntranscend:=Modulet(}, sin[x_]:=Sintx]; cos[x_]:=Cos[x]; 
tan[x_];*Tan[x]; cot[x_]:*Cot[x]; csc[x_]:=Csc[x]; 
sec[x_]:=Sec[x]; exp[x_];=Exptx] ] 

nSo1ve[{x==a_,y==d_:0+c_:1* Log[b_:1*x_])]; = 

{a,N[d+c*Log[b*a]]) 

nSolve[{y==d_:0+c_:l*Log[b_:l*x_],x==a_)]:* 

{a,N[d+c*Log[b*a]]) 

nSolve[{y==dl_:0+cl_:l*Log[bl_:l*x], 
y==d2 :0+c2_:l*Log[b2_:l*x]}]:= 

N[{((bl-cl)/(b2-c2)Exp[dl-d2])-(l/(c2-cl)),dl+cl*Log[bl*x]/. 

x->((bl-cl)/(b2^c2)Exp[dl-d2])^(1/(c2-cl))}] 
nSolvetF_List];={x,y) /.NSolve[F,{x,y}] 






II. FILE CARTPLOT.KA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 

cartplot:=Module[{),allyfunc={);allxfunc={);xlines={); 

inyfunc*() ;endpoints={) ;intersect={) ; 

cltranscend; 

f[{x _ )] :=x; 

type2D*l; 

Whi1e[type2D==1, 

yfunc=Input["Enter the curves, one at a time. When 
finished type n. 

y = ’•]; 

yfunc=yfunc /. sqrt->Sqrt; 
yfunc=yfunc /. pi->Pi; 

If[yfunc=!=n, { 

AppendTo[allyfunc, yfunc], 

AppendTo[myfunc, ToExpression[StringJoin["y==","yfunc"]]]), 
{type2D=0}] 

]; 

While[type2D==0, 

xfunc=Input["Enter lines perpendicular to the x-axis. 

When finished type n. 


X = " ]; 

xfunc=xfunc /. sqrt->Sqrt; 
xfunc=xfunc /. pi->Pi; 

If[xfunc=!=n, { 

AppendTo[allxfunc, xfunc], 

AppendTo[myfunc, ToExpression[StringJoin["x==","xfunc"]]]}, 
{type2D=-l}] 

]; 

xlower=Input["Enter lower limit for x-range of plot."]; 
xlower=xlower /. sqrt->Sqrt; 
xlower=xlower /. pi->Pi; 

xupper=Input["Enter upper limit for x-range of plot."]; 
xupper=xupper /. sqrt->Sqrt; 
xupper=xupper /. pi->Pi; 

ylower=Input["Enter lower limit for y-range of plot."]; 
ylower=ylower /. sqrt->Sqrt; 
ylower=ylower /. pi->Pi; 

yupper=Input["Enter upper limit for y-range of plot."]; 
yupper=yupper /. sqrt->Sqrt; 
yupper=yupper /. pi->Pi; 
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I f [Length[alIxfunc]>»!,{ 

For[i=l, i<=Length[allxfunc], {ival=i; ++i; 

1ine-show[Graphics[{ 

Line[{{allxfunc[[ival]], ylower),{axlxfunct[ival]]/ 

yupper)}] )], AxesLabel->{ "x" , ••y*'}, DisplayFunction 
->Identity]; 

AppendTo[xlines,1ine] )]; 

}] ; 

rmtranscend; 

theplot=Plot[Evaluate[allyfunc],(x, xlower, xupper). 
AxesLabel~>("x","y "), PlotRange“>((xlower,xupper), 

{ylower,yupper)},DisplayFunction->ldentity); 

Show[xlines, theplot,DisplayFunction“>$DisplayFunction]; 

X*• f y= •f 

last=Length[rayfunc]; 

For[ii*l, ii<= last, 

(ival=ii; ++ii; 

For[jj=l, jj<=last, 

{jval=jj; ++jj; 


cltranscend; 

If[jval!=ival && 
(myfunc[[ival,l]]*==y 


inyfunc[[jval,l]]===y), 


ln[x_]:=Log[x]; 

If [inyfunc[ [ival, 1 ] ]=«=y && inyfunc[ [ jval,l] ]=*=y. 

If [ {inyfunc[ [ival] ], Biyfunc[ [ jval] ] } = ! = 
({myfunc[[ival]], inyfunc[ [ jval] ] /. Log[x_]->x)), 
{theroots-{nSolve[ (myfunc[ [ival] ] ,inyfunc[ [ jval] ]} ]} 

} r 


transcend; 

theroots=NSolve [ {inyfunc[ [ival] ], inyfunc[ [jval] ]), {x,y}]; 
theroots=N[{x,y) /. theroots]; 
theroots=DeleteCases[DeleteCases[DeleteCases[ 
theroots, {x_Coinplex,y_) ], {x_,y_Coniplex) ], 

{x_Complex, y_Coinplex} ]; 

If [ Length [ theroots ] "0, { 

theroots=Solve[{myfunc[[ival]], ayfunc[[jval]]},{x,y}]; 
theroots=N[{x,y} /. theroots]; 
theroots=DeleteCases[DeleteCases[DeleteCases[ 
theroots, {x_Coinplex,y_) ], {x_,y_Coinplex} ], 

(x_Coinplex, y Complex) ] 

}]; 


cltranscend; 


If[theroots=!={x,y} && theroots !={) && 

{inyfunc[ [ival] ] ,inyfunc[ [ jval] ] )=== 

(({inyfunc[ [ival] ] ,niyfunc[ [ jval] ])) /. Sqrt[x_]->x) , 

{ 

rmtranscend; 

theroots=Table[FindRoot[myfunc[[ival,2]]- 
inyfunc[ [jval, 2] ]"0, (x, theroots [ [ffl, 1] ]} ], 

{m,Length[theroots])]; 
theroots=x /. theroots; 

theroots=Table[{theroots[[ffl]],(fflyfunc[[ival,2]] /. 
x-> theroots[[ffl]])},{ffl,Length[theroots])]; 
temproots= 

Table[(TrueQ[(Abs[fflyfunc[[ival,2]]-fflyfunc[[jval,2]]] /. 

X -> therootsi[s,1]]) <=.5]),{s,Length[theroots])]; 
theroots==Transpose[ (theroots,temproots) ]; 

theroots=Transpose[DeleteCases[theroots, {_, (False)) ,1] ] [ [1] ]; 
If[theroots*= {)[[!]],theroots={) ]; 

x=. 

)]; 

)]; 

)]; 

rmtranscend; 

If [myfunc[ [ival, 1] ]*=*=y && myfunc[ [ jval, 1] ]===x, 
theroots*N[{(myfunc[[jval,2]],(myfunc[[ival,2]]/. 
x->myfunc[[jval,2]])}}]]; 

If[fflyfunc[[ival,1]]=*=x && myfunc[[jval,1]]*==y, 
theroots=N[((fflyfunc[[ival,2]],(fflyfunc[[jval,2]]/. 
x->myfunc[[ival,2]]))}]]; 

theroots=DeleteCases[DeleteCases[DeleteCases[ 
theroots, {x_Complex,y_}],(x_,y_Complex>], 

{x_Complex,y_Complex}] ; 

theroots=ReplaceAll[theroots, 0. -> 0]; 

If[theroots =!= {x,y} && theroots =!= {),( 
endpoints=Join[endpoints, theroots] )]; 

}]; 

If[Length[endpoints]>0,{ 

endpoints=Select[endpoints,(#[[1]]>=xlower )&]; 
endpoints=Select[endpoints,(#[[!]]<=xupper)&]; 
endpoints=Select[endpoints,(#[[2]]>=ylower)&]; 
endpoints=Select[endpoints,(#[[2]]<=yupper)&]; 

)]; 
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}] ; 

)] ; 

intersect=Union[N[Round[endpoints*1000]/1000],{}]; 


If[Length[intersect]==1, 

Print["The curves intersect at the point: ",f[intersect]]]; 
If[Length[intersect]>=2,{ 

Print["The curves intersect at the points:"]; 

Print[intersect] )]; 
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III. FILE EXAMPLES.HA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 


gr[1]:=Show[shades, 

Plot[ .4 Sin[x] -♦■.2 , {X, .2,1.25), 

OisplayFunction -> Identity, 

AxesLabel -> (’’x", "y"}, Ticks-> None, AspectRatio -> 1, 
PlotRange -> {(0,1.35},{07,.6))], 

Graphics[( 

Text["a", {.2,-.05)], 

TextC’b**, {1.25,-.05)], 

Text["y = f(x)", {.75,.6)] 

)], 

OisplayFunction -> $DisplayFunction]; 

example[1]:=Show [ 

Plot[x^2, {x,0,l}, 

OisplayFunction -> Identity, 

AxesLabel -> {"x","y”)» 

Ticks -> {{1},{1)}, 

AspectRatio -> 1, 

PlotRange -> {{-.06,1.1),{-.1,1))], 

Graphics[{ 

Line [{{1,0), {1,1}}], 

Text["0”,{.04,-.08)], 

Text ["y » x", {.5,.8)], 

Text [”2",{.615,.85)] 

)], 

OisplayFunction -> $DisplayFunction]; 

example[2];=Show[ 

Plot[ Cos[x], {X,-Pi/2, Pi/2), 

OisplayFunction -> Identity, 

AxesLabel -> {"x",”y")f 
AspectRatio->.5, 

Ticks -> None, 

PlotRange -> { {-1.65,1,65), {-.3, 1.2)}], 

Graphics[{ 

Text[ " y = cos x", {1.4,1.05)], 

Text[ '•-Pi/2", {-Pi/2,-.2)], 

Text[ "Pi/2", {Pi/2,-.2)] 

)], 

OisplayFunction -> $DisplayFunction]; 
gr[2]:=Show[ 

Plot[.5*{x-1)^3 +.2 ,{x,.5,1.5), OisplayFunction -> Identity, 
AxesLabel -> {"x","y"), Ticks-> None, AspectRatio ->1, 
PlotRange -> {{0,2),{-.085,.265))], 

Graphics[{ 
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Line[{ (.5, 0}, {.5, (.5-1)^3+.325) )], 

Line[{ {1.5, 0), (1.5, .27))], 

Line[( (1.68,0), (1.68,.088))], 

Line[( (1.68,.13), (1.68,.2))], 

Line[( (.5,.2},(1.5,.2))], 

Line[( (1.58,.2), (1.78,.2))], 

Line[( (.5, -.08),(.5,-.0525))], 

Line[( (1.5,-.08), (1.5,-.0525))], 

PointSize[.01], Table[Point[(l,n .005)],(n, 0, 39)], 

Text fa'*, (.5, -.025)], 

TextCb", (1.5, -.025)], 

Textt"c", (1, -.025)], 

Textt" f(c)", (1.68, .11)], 

Textt''0'', (-.075,-.025) ], 

Text["<- (b - a) ->«, (1.005, -.0675)], 

Textfy * f (X)", (1, .25)], 

Textt” V "# (1.68,. 01 )], 

Text(" /\ (1.68, .19)] 

DisplayFunction -> $DisplayFunction]; 

example[3]:=Show[Plot[(Sqrt[4-x^2],-Sqrt[4-x^2]),(x,-2,2), 
DisplayFunction ->Identity, 

Ticks -> ((-2,0,2),()), 

AxesLabel -> ('’x”, "y"), 

AspectRatio -> Automatic, 

PlotRange -> ((-2.15,2.2),(0,2.1))], 

Graphics[{ 

Text[”y * 4 - X ", (2,2.35)], 

Text["2",(2.8,2.65)] 

)]/ 

DisplayFunction -> $DisplayFunction]; 

example[4]:=show[Plot[2x,(x,l,4), 

DisplayFunction ->Identity, 

Ticks -> ((1,4),()), 

AxesLabel -> ("x'’,"y"), 

AspectRatio -> l, 

PlotRange -> {(0,4.6),(0,8))], 

Graphics[{ 

Text["y = 2x", (2,7)], 

Text[»(l,2)",(.5,2)], 

Text["(4,8)",(4.6,8)] , 

Line[({l,0),{l,2))], 

Line[{(4,0),(4,8))] 

DisplayFunction -> $DisplayFunction]; 

example[5];=Show[Plot[(Cos[x], -Sin[x]),(x,-.1,1.6), 
DisplayFunction -> Identity, AxesLabel -> ("x","y"). 

Ticks -> ((0),(-1,0,1)), AspectRatio -> Automatic, 
PlotRange -> ((-.15,1.8),(-1.2,1.1))], 
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Graphics[( 

Line[({Pi/2,0),{Pi/2,-1))], 

Text["Pi/2", {1.75,.2)), 

Texti^y = cos x", {1,1.1)], 

Texti^y « -sin x", {.75,-1.1)] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[6]:=Show[shadel,Plot[l-x^(-2),{x,.5,5), 
DisplayFunction -> Identity, 

AxesLabel -> {”x'',"y"). Ticks -> {{1,4),{!)), 

AspectRatio -> Automatic, 

PlotRange -> {{-.1,5),{-1.75,1.5))], 

Graphics[{ 

Line[{{0,l),{5.5,l))], 

Line[{{l,0),{l,l))], 

Line[{{4,15/16),{4,1))], 

Text["y = 1",{2,1.3)], 

Textt”y = 1 - X",{1.6,-!)], 

Textt"-2", {2.6,-.75)] 

}], 

DisplayFunction -> $DisplayFunction]; 

example[7];=Show[Plot[{-x, 2-x^2),{x,-l,2), 
DisplayFunction -> Identity, AxesLabel -> {"x","y”)» 
Ticks -> {{-1,0,1,2),{-2,0,2)), AspectRatio -> Automatic 
PlotRange -> {{-1.95,2.2),{-2.4,2.1))], 

Graphics[{ 

Text["(-l,l)", {-1.65,1)], 

Text["(2,-2)", {2,-2.2)], 

Text[”2'', {2.35,1.9)], 

Text[''y = 2 - X ”, {1.7,1.7)], 

Text["y = - X", (.9,-1.8)] 

) 1 , 

DisplayFunction -> $DisplayFunction]; 

example[8]:=Show[Plot[{x^3+l,x^2+x),{x,-1.5,1.2), 
DisplayFunction ->Identity, AxesLabel -> {”x”,”y”). 

Ticks -> {{1),{2)), AspectRatio -> Automatic, 

PlotRange -> {{-1.75,1.5),{-.5,2.3))], 

Graphics[{ 

Text["y = X + 1”,{-.7,1.2)], 

Text["3",{-.625,1.35)], 

Text[”2”,{1.175,.65)], 

Text[”y = X + x",{l.l,.5)], 

Text[”(-l,0)”,{-1.5,-.2)], 

Text[”(l,2)”,{1.4,2)] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[9]:=Show[Plot[{4+x^2,12-x^2),{x,-2.2,2.2), 
DisplayFunction -> Identity, Ticks -> None, 








AxesLabel -> , AspectRatio -> Automatic, 

PlotRange -> {{-4.5,4.5), {2,13})], 

Graphics[( 

Textt**(-2,8) •',(-3.2,8}], 

Text[''(2,8)",{3.1,8)], 

TextCy = 12 - X", (-3,12.5)], 

Text[••2'’, (-1.125,12.95)], 

TextCy = 4 + x", (3,3.5)], 

Text[''2", (4.725,3.95)] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[10]:=Show[Plot[Sqrt[x],{x,0,4), 

DisplayFunction -> Identity, 

AxesLabel -> ("x","y"), Ticks -> ((2,4),(1,2)), 
AspectRatio -> .75, 

PlotRange -> ((-.2,4.5),(-.2,2.2))], 

Graphics[( 

Line[((2,0),(2,Sqrt[2])}], 

Line[((2,0),(4,2))], 

Text["2",(1.7,2.2)], 

TextL"y = X", (1.85,2.1)], 

Text["(4,2)",(4.4,2)], 

Text["y = 0", {1,-.15)], 

Text["y = X - 2", (4,.9)] 

}], 

DisplayFunction -> $DisplayFunction]; 

example[11]:=Show[shade2, Plot[2-x^2,(x,0,Sqrt[2]), 
DisplayFunction -> Identity, 

AxesLabel -> {"x","y"), 

Ticks -> ((1),(1,2)), AspectRatio -> Automatic, 
PlotRange -> ((0,2),(-.07,2.2))], 

Graphics[( 

Text[’'y = 2 - x", (.65,2.1)], 

Text["2",(1,2.2)], 

Text["y = X", (1.65,1.25)], 

Text["y = 0",(.55,-.!)], 

Line[({0,0),(1.7,1.7))], 

Line[{(2.1,2.1),(2.3,1.55)}], 

Line[((l,0),(1,1))] 

}], 

DisplayFunction -> $DisplayFunction]; 

example[12]:=Show[shade!,Plot[(Sqrt[x],-Sqrt[x],6-x), 
(x,-.l,ll}, 

DisplayFunction -> Identity, AspectRatio -> Automatic, 
AxesLabel -> ("x'*,"y"}. Ticks -> None, 

PlotRange -> ((-.1,11),(-4.2,4.2))], 

Graphics[( 

Line[((0,0),(9,0))], 

Text["(9,-3)",{10,-2.2)], 
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Text[«(4,2)«,{4.2,3.2)], 

Text["y « X", {2.2,-3)] 

Text["2",{1.85,-2.67}], 

Text["y * 6 - x",{7.2,1.!)] 

}]/ 

OisplayFunction -> $DisplayFunction]; 

gr[3]:»Show[Plot[4(x+.l)^3,{X,.2,.5), 

OisplayFunction -> Identity, 

AspectRatio -> .4, PlotRange -> {{-.1,1.2},{0,1}}, 
AxesOrigin -> {-.1,0}, 

AxesLabel -> {"x","y"},Ticks -> None], 

Plot[4(X-.4) '3,{X,.7,1}, OisplayFunction -> Identity, 
AspectRatio ->.5, PlotRange -> {{0,1},{-.1,1.3}}], 

Graphics[{ 

Line[{{.485,.81},{.99,.81}}], 

Line[{{.248,.175},{.756,.175}}], 

TextCa", {-.15, .175}], 

Text["b",{-.15,.81}], 

Text["x = g(y)-,{.2,.95}], 

Text["x = f(y}",{1.1,.2}] 

}], OisplayFunction -> $OisplayFunction]; 

gr[4]:=Show[Plot[2.1-(x+.3)^2,{X,.59,1.1}, 

OisplayFunction -> Identity, 

AspectRatio ->.4,AxesLabel-> {"x","y”},Ticks -> None, 
AxesOrigin -> {-.1,0}, 

PlotRange ->{{-.1,1.3},{0,1.4}}], 

Plot[ (X+.5)''3-.2, {X, .2, .66}, OisplayFunction -> Identity, 
Ticks -> None, 

AspectRatio ->.4, PlotRange -> {{0,1.3},{-.1,1.3}}], 
Graphics[{Line[{{.25,.2},{1.08,.2})], 

Text[''a", {-.15, .2}], 

Text["b",{-.15,1.235}], 

Text["x = f(y)",{l.2,.9}], Text["x = g(y)”,{.2,1}]}], 
OisplayFunction -> $0isplayFunction]; 

gr[5]:=Show[Plot[{-Sqrt[x^2-.5]+1.5,Sqrt[x^2-.5]+1.5},{X,. 
711,1.5), 

OisplayFunction -> Identity, AspectRatio ->.5, 

AxesLabel -> {”x”,”y”}. Ticks -> None,AxesOrigin -> {-.1,0}, 
PlotRange -> {{-.1.2.55},{0,3}}], 

Plot[{-Sqrt[3-x^2 3 ^1.5,Sqrt[3-x^2]+1.5},{x,1.2,1.732051), 
OisplayFunction Identity, AspectRatio -> .5, 

Ticks -> None],Graphics[{ Text[”x = f(y)”,{2.2,.5}], 

Line[{{1.73205,1.45},{1.73205,1.55})], 

Text["a",{-.25,.38)], Text[”b”,{-.25,2.62}], 

Text["x = g(y)”,{.55,2.7}]}], 

OisplayFunction -> $OisplayFunction]; 
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example[13]:=Show[Plot[Sqrt[x],{x,0,4}, 

DisplayFunction -> Identity, 

AxesLabel -> {"x","y"), Ticks -> {{2,4),{1,2}), 

AspectRatio -> Automatic, 

PlotRange -> {{-.2,4.5),(-.45,2.2)}], 

Graphics[{ 

Linet{{2,0),{4,2))], 

Text[''2", {2.35, 2.3)], 

Text["y = X", {1.8,2.!)], 

Text[”(4,2)«,(4.4,2)], 

Text[”y = 0", {l,-.3)], 

Text["x = y + 2", (4,.9)] 

)]/ 

DisplayFunction -> $DisplayFunction]; 

example[14]:=Show[shade2, Plot[2-x^2,{x,0,Sgrt[2]), 
DisplayFunction -> Identity, 

AxesLabel -> {’•x”, "y"), 

Ticks -> {{1),{1)), AspectRatio -> Automatic, 

PlotRange -> {(0,2.25),(-.22,2.05))], 

Graphics[{ 

Text[''(l,l)”,(1.3,1)], 

Text[''2", (1.5, .67)] , 

Text["x = 2 - y”, (1.7,.6)], 

Text[«y = X", (.4,.8)], 

TextC^y = 0'*, (.55,-.15)], 

Line[((0,0),(1.5,1.5))] 

)]. 

DisplayFunction -> $DisplayFunction]; 

example[15]:=Show[shade!,Plot[(Sqrt[x],-Sqrt[x],6-x),(x,-.l, 
11), 

DisplayFunction -> Identity, AspectRatio -> Automatic, 
AxesLabel -> ("x","y”), Ticks -> None, 

PlotRange -> ((-.1,11),(-4.2,4.2))], 

Graphics[{ 

Line[((0,0),(9,0))], 

Text[''(9,-3)'’, (10,-2.2)] , 

Text['’(4,2)", (4.2,3.2)] , 

Text["x = y",(2.8,-3.2)], 

Text["2", (3.875,-2.8)], 

Text[''x = 6 - y",(7.5,.9)] 

> 3 , 

DisplayFunction -> $DisplayFunction]; 

gr[6];=Show[shadell,Plot[Sqrt[l-x^2],(x,.59,.975}, 
DisplayFunction -> Identity, 

AxesLabel-> ("x", ’’y"), AspectRatio -> Automatic, 

Ticks -> None,PlotRange -> {{0,1.2),(0,.8))], 

Plot[Sqrt[.5-x^2],(X,.4375,.7071),DisplayFunction -> Identity, 
Ticks -> None,AspectRatio -> 1], 

Plot[Sqrt[.25-x^2],(X,.45,.5),DisplayFunction -> Identity, 
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Ticks -> None,AspectRatio ->1], 

Graphics[{ 

TextCr-fCQ)", {1.15, .255)], 

TextC^b",{.575,.5}], 

Text["a'',{.55,.125)], 

Text[''A",(.9,.55)], 

TextC^B", {.625,.865}], 

Line[((0,0),{l,.5))], 

Line[{{0,0},{.7,.88}}] 

>], 

DisplayFunction -> $DisplayFunction]; 
example[16]:=Show[PolarPlot 

[2(l+cos[Q]),(Q,0,2Pi), DisplayFunction -> Identity, 
AspectRatio -> Automatic, AxesLabel -> {"x","y"), 

Ticks -> None,PlotRange -> {(-1,6},{-2.65,3})], 

Graphics[{ 

Line[((0,0),{2.414, 2.414)}], 

Line[{{4,0},{4.25,.225}}], 

Textt''2", {-.2,2.1}], 

Text[''-2'', {-.25,-2.2}] , 

Text t "0=0,2Pi'' ,{5.1,.35}], 

Text['*4", {4.2,-.25}], 

Text[«P{r,Q)«, {2.6, 2.75}], 

Text["r'', {1.2, 1.45}], 

Text["r= 2{l+cosQ}",{2.5,3.35}] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[17]:=Show[shadel2,PolarPlot[{l+Sin[Q]},{Q,0,2Pi}, 
DisplayFunction -> Identity, Ticks -> None, 

AspectRatio -> Automatic, 

AxesLabel -> {"x'',"y"}, PlotRange -> {{-1.8,2.5}, {-1.1,2.3}} ], 
Graphics[{ 

TextCr = 1 + sinQ", {-.95,2.2}], 

Text["Q=Pi/4", {1.65,1.65}], 

Text[''Q=-Pi/4'', {.85,-.85}], 

Line[{{0,0},{1.4,1.4}}], 

Line[{{0,0},{.7,-.7}}] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[18];=Show[shade7,PolarPlot[{2 Cos[Q]-H),{Q,0,2Pi}, 
DisplayFunction -> Identity, Ticks ->None, 

AspectRatio -> Automatic, AxesLabel -> {”x”,’'y"}» 

PlotRange-> {{-1.9,4},{-2,2.3}}], 

Graphics[{ 

Line[{{-.75,1.3},{.5,-.866}}], 

Line[{{-.75,-1.3},{.5,.866}}], 

Line[{{3,0},{3.2,.3}}], 

Line[{{l,0},{1.2,.3}}], 

Line[{{0,0},{3,0}}], 
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Textt^r » 2 COSQ + 1", {1.75, 2.15)], 

Text[''Q»2Pi/3'', {-1.2,1.45)], 

Text["QMPi/S", {-1.2,-1.45)] , 

Textt'*Q=Pi'*, {1.65, .4)], 

Text["Q»0",{3.55,.4)] 

)]. 

DisplayFunction -> $DisplayFunction]; 

example[19]:BShow[shadelO,PolarPlot[{Sin[3Q]),{Q,0,2Pi), 
DisplayFunction -> Identity, Ticks -> None, 

PlotRange -> {{-.9,1.2),{-1.05,1.1)), 

AspectRatio -> Automatic, AxesLabel -> {"x","y")], 
Plot[Sqrt[l. l-x''2], {X, .908,1.05),DisplayFunction -> Identity, 
Ticks -> None, AspectRatio -> Automatic], 

Graphics[{ 

Line[{{0,0),{1.039,.6))], 

Line[{{0,0),{.6,1.039))], 

Text["r =sin3Q'', {-.6,.9)], 

Text[''Q=Pi/3", {.6,1.15)], 

Text [”Q=Pi/6'’,{ 1.15,. 725)], 

Text["Q=0'', {!,-.!)], 

)}. 

DisplayFunction -> $DisplayFunction]; 

example[20]:=Show[PolarPlot[ {Sqrt[Cos[2Q]]), {Q,0,2Pi), 
DisplayFunction -> Identity, AxesLabel -> {"x'*,”y'*). 

Ticks -> None, AspectRatio -> Automatic, 

PlotRange -> {{-1.1,1.1),{-.8,1.1))], 

Graphics[{ 

Line[{{-.6,-.6),{.6,.6))], 

Line[{{-.6,.6),{.6,-.6))], 

Text[''r = COS2Q", {-.4, .95)], 

Text["2'', {-.6,1.05)1, 

Text["Pi/4",{.7,.7)], 

Text["-Pi/4",{.7,-.7}], 

Text[''3Pi/4'', {-.75, .7)] , 

Text["5Pi/4",{-.75,-.7)], 

TextC'l", {1.05,-.!)] 

}], 

DisplayFunction -> $DisplayFunction]; 

gr[7];=Show[shade6,Plot[Sqrt[l-x^2],{X,.555,.896), 
DisplayFunction -> Identity, AxesLabel -> {"x","y”)» 

Ticks -> None, AspectRatio -> Automatic], 

PIot[.5x,{X,0,1.275),DisplayFunction -> Identity, 

Ticks -> None], 

Plot[1.5x,{x,0,.665),DisplayFunction -> Identity, 

Ticks ->None], 

Plot[Sqrt[.75-x^2+X],{X,.658,1.27), 

DisplayFunction -> Identity, Ticks -> None], 

Graphics[{ 

Text[”rl",{.65,.65)], 
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Text["r2",{l,.95)], 

Text[''Q=*a'*, {1.2, .475}], 

TextC^Q-b",{.45,.95)] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[21]:^Show[shade5,shade4, 

PolarPlot[l,(Q,0,2Pi},AxesLabel -> {"x","y"), 

Ticks -> None, AspectRatio -> Automatic, 

DisplayFunction -> Identity], 

PolarPlot[.85,(Q,0,1.3), 

Ticks -> None, AspectRatio -> Automatic, 

DisplayFunction -> Identity], 

PolarPlot[(l-Cos[Q]),{Q,0,2Pi), DisplayFunction -> Identity, 
Ticks -> None, AspectRatio -> Automatic], 

Graphics[( 

Line[((0,-1},(.2,-1.3))], 

Line[{{0,l},{.2,1.3}}], 

Line[((0,0),{.2,.98}}], 

Line[{{0,0},(.3,.954})], 

Line[{{0,0),{l,0})], 

Text [ •»Q=-Pi/2" , { . 5, -1.4 } ] , 

Text[«Q=Pi/2",{.4,1.4}], 

Text[”Q«,(.6,.4}], 

Text['’r2 = 1", (.55,1.15}], 

Textt''rl = 1 - cosQ", {-1.1,1.5}] 

)]/ 

DisplayFunction -> $DisplayFunction]; 

example[22] :=Show[shade9,PolarPlot[ {3CostQ] , 
l+Cos[Q]},{Q,0,2Pi}, 

DisplayFunction -> Identity, AspectRatio -> Automatic, 
AxesLeUsel -> {'*x”,"y"}. Ticks -> None, 

PlotRange -> {{-.3,3.75},(-2.5,2.5}}], 

Graphics[{ 

Line[{(0,0},{1.25,2.16}}], 

Line[{{0,0},{3,0}}], 

Line[{{0,0},{1.25,-2.16}}], 

Line[{{1.6,-1},{2.2,-1.9}}], 

Line[{{1.9,1.8},{1.7,1.5}}], 

Text["Q=Pi/3",{1.25,2.3}], 

Text["Q=-Pi/3",{1.25,-2.3}], 

Text[''r2 = 3COsQ'', {2.7,1.95}], 

Text["rl = 1 + cosQ",{3.05,-2.05}] 

)], 

DisplayFunction -> $DisplayFunction]; 

example[23];=Show[PolarPlot[{l-Cos[Q]},{Q,0,2Pi}, 
DisplayFunction->Identity, 

Ticks->None, AxesLabel->{"x","y”},AspectRatio->Automatic, 
PlotRange->{{-2.75,1},{-1.3,1.3}}], 

Plot[Sqrt[.l-x^2],{x,-.161,.318},DisplayFunction->Identity, 
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Ticks-> None,AspectRatio->Autoinatic ], 
Graphics[{ 

Linet{{0,0),{-.75,1.3))], 

Text("Q'*,{.35,.2)], 

Text[''r'',(-.55,.55)], 

TextCr - l-cosQ", {-2.2,1.2}], 

Text[••P(r,Q)", {-.85,1.45)3 , 

Text[••-2'*, {-2.15,-.15)3, 

TextCl", {-.1,1.2)3 

) 3 , 

DisplayFunction->$DisplayFunction3; 
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IV. FILE mTBGRAL.MA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 


BeginPackage["integral'"]; 


graph::usage = "Type graph, to invoke a program which 
plots a region, finds the limits of integration, and 
evaluates the integral." 


SetOptions[Plot, AxesLabel -> {"x","y"}]; 
SetOptions[ParametricPlot, PlotPoints -> 600]; 
SetOptions[Plot,, PlotPoints->150] ; 

Of f[Drop::drop]; 

Off[EndPackage::noctx] 

Off[FindRoot::cvnwt]; 

Off[General::spell]; 

Off[General::spelll]; 

Off[General::stop]; 

Off[Get::noopen]; 

Off[Graphics::gprim]; 

Off[Graphics::gptn]; 

Off[Greater::nord]; 

Of f[Infinity::indet]; 

Off[NIntegrate::ncvb]; 

Off[NIntegrate::nlim]; 

Off[NIntegrate::slwcon]; 

Off[ParametricPlot::pptr]; 

Off[Part::partd]; 

Off[Part::partw]; 

Off[Plot::plln]; 

Off[Plot::plnr]; 

Off[Power::infy]; 

Of f[Remove::remal]; 

Off[ReplaceAll::rmix]; 

Off[ReplaceAll::reps]; 

Off[Show::gtype]; 

Off[Solve::ifun]; 

Off[Solve::tdep]; 

Off[Unset::norep]; 


Needs["Graphics'Master'"]; 
«c:polarl.ma; 

«c:polar2 .ma; 

«c: cartplot. ma; 

«c: inverse. ma; 
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«c: examples. ma; 

«c; shades, ma; 

«c:cart.ma; 

«c: polplotl. ma; 

«c: polpl ot2. ma; 

Begin["'Private'"]; 

graph:* 

Module[{}, 

labexamplesInput[''To display a graph that is in the lab 
book type in the example number, otherwise type n."]; 

If[labexample =!= n, example[labexample], { 

type2D=Input["Type 1 for cartesian coordinate system, or 
type 2 for polar coordinate system."]; 

If[type2D==l, { 

plotonly=Input["To just plot curves type y, otherwise 
type n."]; 

If[plotonly===y, cartplot, cart] 

)]? 


If[type2D==2, { 

plotonly=Input["To just plot curves type y, otherwise 
type n."]; 

If[plotonly===y, { 

typepolplot=Input["Type 1 for plot of area between the 
origin and a polar curve, or type 2 for plot of area 
between two polar curves."]; 

If[typepolplot==l,polplotl]; 

If[typepolplot==2,polplot2]}, { 

typepolar=Input["Type 1 for area between the origin and 
a polar curve, or type 2 for area between two polar curves."]; 
If[typepolar==l, polarl]; 

If[typepolar==2, polar2]; 

)]; 

)]; 

)]; 

] 

End[]; 

EndPackage[]; 


127 





▼ . FZLB mVSRSB.MA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 

inverse:^Module[{}, 
ln[x_]z-Logtx]; 

If [myfunc=“(inyfunc /. Log[x_]->x), { 

transcend; 
f[(x_}];» x; 

allycoord=Transpose[Join[f[thelines]]][[2]]; 
ycoord=sUnion[ycoord, {Max[all vcoord]), {Mintallycoord]) ]; 
ycoord=N[Union[Round[ycoord*100]]/10oi; 
ycoord=Union[ReplaceAll[ycoord, 0.->0],()]; 
xfuncs-Transpose[PartitionJJoin[myfunc, thelines ], 

Length[myfunc]]]; 

yfunc=Complement[xfuncs,Cases[xfuncs, {y==b_Integer,c_) } 
(y==b_Real,c_)]]; 

For[ww=l, w<=Length[yfunc], {wval=ww; ++ww; 

If[yfunc[[wval,l,l]]===y, 

yfunc[ [wval, 1] ]=ToExpression[StringJoin[''x==**, 

"(X/. Solve[yfunc[[wval,1,2]3==y, x, 

InverseFunctions-> True])[[1]]"]] ]; 

If[yfuncC[wval,l]]==*True || 

(yfunc[[wval,1,2]] /. (y_)^(1/3)->y)=!«yfunc[[wval,1,2]], 

ww=99] 

)]; 

If[WW<99, { 

revlines=Transpose[yfunc][[2]]; 

For[vv=l, vv<=I,ength[ycoord]-l, {vval=vv; ++vv; 
yy=(ycoord[ [wal] ]+ycoord[ [wal+l] ] )/2; 

integrfunc=Select[revlines,((yy> Min[Transpose[#] [[2]]]) && 
(yy< Max[Transpose[#] [[2]]]))&]; 

If[Length[integrfunc]>2,(vv=99; ww=99)]; 

)]; 

)]; 

If[ww<99, { 

Pause[8]; 

yintegr=Input["Type y if you wish to integrate the same 
problem with respect to y, otherwise type n.”]; 

If[yintegr===y,( 

Print[" "]; 

Print[••**********************************************•• j; 

invfunc=Transpose[yfunc][[1]]; 
tempinvfunc=invfunc; 
revlines=Transpose[yfunc][[2]]; 


128 







For[zz*l, zz<=Length[ycoord]-l, (zval=zz; ++zz; 
yy=(ycoord[[zval]]+ycoord[[zval+1]p/2; 

integrfunc=Select[revline8,((yy> Min[Transpose[#] [[2]]]) && 
{yy< Max[Transpose[#] [[2]i]))6]; 

invf unc=teinpinvf unc; 

If[Length[integrfunc]*=2,{ 
cl=Drop[Select[integrfunc[[1]] , 

(#[[2]]>=ycoord[[zval]] && #[[2]]<=ycoord[[zval+l]l)6],l]; 
c2*Drop[Select[integrfunc[[2]] , 

(#[[2]]>=ycoord[[zval]] && #[[2]]<=ycoord[[zval+l]])S],-l]; 
avgcurvel=(Transpose[cl][[1,1]]+Transpose[cl][[1,-1]])/2; 
avgcurve2=(Transpose[c2][[1,1]]+Transpose[c2][[1,-1]])/2; 

adjcl=Transpose[Reverse[Transpose[cl]]]; 
inaxcl=Max[Union[{adjcl[ [-1,1] ]), {adjcl[ [1,1] ]) ] ] ; 
mincl=Min[Union[{adjcl[[-l,l]]},(adjcl[[l,l]])]]; 
curvel=NIntegrate[Fit[adjcl,{l,y,y^2),y], 

{y,inincl,inaxcl) ]/(maxcl-mincl) ; 
adjc2=Transpose[Reverse[Transpose[c2]]]; 
inaxc2=Max[Union[ {adjc2[ [-1,1]]), {adjc2[ [1,1]]}]]; 
ininc2=Min [Union [ {ad j c2[[-l,l]]),( ad j c2 [ [ 1,1 ] ]} ] ]; 
curve2=NIntegratetFit[adjc2,{l,y,y^2),y], 

(y ,minc2 ,maxc2) ]/(inaxc2-ininc2) ; 

}]; 

If[Length[integrfunc] -~ 2 ,{ 

If[avgcurve2>avgcurvel,{integrfunc=Reverse[integrfunc]; 
tempcurve=curve2; curve2=curvel; curve l=teinpcurve; 
teinpcurve=avgcurve2; avgcurve2=avgcurvel; avgcurvel=teinpcurve) ] 
)]; 

If[Length[integrfunc]==2,{ 

placel=Position[revlines, integrf^nc[[1]]][[1]][[1]]; 
place2=Position[revlines, integrfunc[[2]]]["i]][[1]]), 
{theplace=Position[revlines, integrfunc[[l]]][[l]][[l]])]; 

If[Length[integrfunc]==2,{ 

If[Replace[invfunc[[placel,2]],Sqrt[x_]->0] =1= 
invfunc[[placel,2]], 

I f [ curveKavgcurvel, 

invfunc[ [placel,2] ]=negbranch[invfunc[ [pl=icel,2] ] ] ] 

] 

)]; 

If[Length[integrfunc]==2,{ 

If[Replace[invfunc[[place2,2]],Sqrt[x_]->0] =!= 
invfunc[[place2,2]], 

If[curve2 <avgcurve2, 

invfunc[[place2,2]]=negbranch[invfunc[[place2,2]] ] ] 

] 

)] ? 
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y». ; 

If[Length[integrfunc]==2, { 

area=NIntegrate[invfunc[[placel,2]]-invfunc[[place2,2]], 

{y, ycoord[[zval]],ycoord[[zval+1]])]},{ 

area=NIntegrate[Evaluate[2bothbranches[invfunc[[theplace,2]] 

]], 

(y I ycoord[[zval ]], ycoord[[zval+l]])])]; 

If[Length[integrfunc]==l, { 
right=invfunc[[theplace,2]]; 
left=negbranch[invfunc[[theplace,2]] ]), 

{right=invfunc[[placel,2]];left=invfunc[[place2,2]])]; 

Print[” ”]; 

If[zval==l &6 Length[ycoord]==2, { 

Print["The limits of integration are y = 
ycoord[[zval]]," and y = ", ycoord[[zval+1]] ]; 

Print["The right curve is x = ", right]; 

Print["The left curve is x = ", left]; 

Print["The area of the region is ", Abs[area] ] }]; 

If[zval==l && Length[ycoord]>=3, { 

Print["The limits of integration for the first region"]; 
Print["are y = ",ycoord[[zval]], 

" and y = ",ycoord[[zval+1]]]; 

Print["The right curve for the first region is x » ", right]; 
Print["The left curve for the first region is x = ", left]; 
Print["The area of the first region is ", Abs[area] ]; 
areal=Abs[area] )]; 

If[zval==2 && Length[ycoord]>*3, { 

Print["The limits of integration for the second region"]; 
Print["are y = ",ycoord[[zval]], 

" and y = ",ycoord[[zval+1]]]; 

Print["The right curve for the second region is x = ", right]; 
Print["The left curve for the second region is x = ", left]; 
Print["The area of the second region is ", Abs[area] ]; 

If[Length[ycoord]==3, 

Print["The area of the total region is ", areal+Abs[area] ], 
area2=Abs[area] ] 

)]? 

If[zval==3 && Length[ycoord]>=4, { 

Print["The limits of integration for the third region"]; 
Print["are y = ",ycoord[[zval]], 

" and y = ",ycoord[[zval+1]]]; 

Print["The right curve for the third region is x = ", right]; 
Print["The left curve for the third region is x = ", left]; 
Print["The area of the third region is ", Abs[area] ]; 

If[Length[ycoord]==4, 

Print["The area of the total region is ", areal+area2+ 

Abs[area]], area3=Abs[area] ] 
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If[zval*=4 && Length[ycoord]*=5, { 

Print["The limits of integration for the fourth region"]; 
Print["are y ® ", ycoord[[zval]], 

" and y * ",ycoord[[zval+1]]]; 

Print["The right curve for the fourth region is x = ", right]; 
Print["The left curve for the fourth region is x = ", left]; 
Print["The area of the fourth region is ", Abs[area] ]; 

Print["The area of the total region is ", areal+area2+area3+ 
Abs[area]] 

}]; 

)]; 

}]; 

}]; 

)]; 

] 

gl[(b_ + Sqrt[a_])*c_]:=(b - Sqrt[a])*c; 
g2[(b_ + d_*Sqrt[a_])*c_]:=(b -d Sqrt[a])*c; 
g3[ b_ + d_*Sqrt[a_] ];= b -d Sqrtfa]; 
g4[ Sqrt[a_] *c_];= -Sqrtfa] *c; 

g5[ Sqrt[a_] ]:= -Sqrtfa]; 

gef b_ + Sqrtfa_] ]:= b -Sqrtfa]; 

hl[(b_+ Sqrt[a_]) *c_] := Sqr1:[a]*c; 

h2f(b_ + d_*Sqrtfa_])*c_]:*d*Sqrtfa]*c; 
h3f b_ + d_*Sqrtfa_] ]:*d*Sqrtfa]; 
h4f Sqrtfa_] *c_];= Sqrtfa]*c; 

hsf Sqrtfa_] ];= Sqrtfa]; 

hef b_ + Sqrtfa_] ]:= Sqrt[a]; 

negbranch[x_]:= 

If[Cases[{X},(b_ + Sqrt[a_])*c_]=!={),gl[x], 

Iffcasesf{X},(b_ + d_*Sqrtfa_])*c_]=!={},g2fx], 

Iffcasesf(x), b_ + d_*Sqrtfa_] ]=!={>,g3fx], 

If[Cases[{x), Sqrt[a_] *c_]=l={},g4[x], 

If[Cases[{X}, Sqrt[a_] ]=!={),g5[x], 

If[Cases[{x), b_ + Sqrt[a_] 3=i={},g6[x]]]]]]]; 

bothbranches[x_];= 

If[Cases[{x),(b_ + Sqrt[a_])*c_]=!={},hl[x], 

Iffcasesf{X},(b_ + d_*Sqrtfa_])*c_]=!={),h2fx], 

Iffcasesf{X}, b_ + d_*Sqrtfa_] ]=l=(},h3fx], 

Iffcasesf{X}, Sqrtfa_] *c_]=!={),h4fx], 

Iffcasesf{X}, Sqrtfa_] ]=!=(},h5fx], 

If[Cases[{x), b_ + Sqrt[a_] ]=!={),h6[x]]]]]]]; 













VI. riLB POLARl.MA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 

polarl:= 

Module[{ 

func, limitskno%ni, theroots, lowerlimit, 
upperlimit, newcurve, area, totalarea*{}, regions={}, 
curvesknown, fund, func2, temp 

}, 

Q=. ; 

cltranscend; 

func=lnput["Enter the polar curve which bounds the region. 


r = " ] ; 

func=func /. sqrt->Sqrt; 
func=func /. pi->Pi; 
func=func /. q->Q; 

revfunc=Replace[func, Sqrt[x_] -> x]; 
limitsknown=Input["If the limits are known type y, 
otherwise type n."]; 

If[1imitsknown===y, 

{lowerlimit=Input["Input lower limit."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit^lowerlimit /. pi->Pi; 

While[lowerlimit <-N[Pi], {lowerlimit=Input["The lowerlimit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

}]; 

upperlimit=Input["Input upper limit."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

)]; 

If[N[upper1imit]<N[lower1imit],lower1imit=N[lower1imit-2 Pi]] 
{lower1imit=0.; upperlimit= 2Pi)]; 
rmtranscend; 

max=Max[DeleteCases[Table[func /. Q -> values, 

{values, 0, 2Pi, .25)],x_Complex]]; 

rayl=N[(func /. Q->lowerlimit)]; 
ray2=N[(func /. Q->upperlimit)]; 
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llequivalent«N[Mod[lowerliinit, 2Pi]]; 
ulequivalentsN[Mcxi[upperliinit, 2Pi] ]; 

If[lowerlimit < 0, llec[uivalent= 2 N[Pi]+lowerlimit, 
llequivalent* lowerlimit]; 

If[upperlimit < 0, ulequivalent» 2 N[Pi]+upperlimit, 
ulequivalent« upperlimit]; 

Q=. ; 

If[(rayl*=ray2) && (llequivalent=*ulequivalent)|| 

(rayl==0 && ray2’==0), { 

PolarPlot[func,{Q, topi[lowerlimit], topi[upperlimit]), 
AxesLabel -> ("x*',"y''}, AspectRatio -> Automatic, 

PlotRange-> {{-1.02 max, 1.02 max),{-1.02 max, 1.02 max}}]}, 
{ 

Show[PolarPlot[func,{Q, topi[lowerlimit], topi[upperlimit]}, 
DisplayFunction -> Identity, AxesLabel -> {"x”,"y"}, 
AspectRatio -> Automatic, PlotRange-> {{-1.02 max,1.02 max}, 
{-1.02 max, 1.02 max}}], 

Graphics[{ 

Line[{{0,0},{Cos[lowerlimit]*rayl, Sin[lowerlimit]*rayl}}], 
Line[{{0,0},{Cos[upperlimit]*ray2, Sin[upperlimit]*ray2}}], 

)]. 

DisplayFunction-> $DisplayFunction]}]; 
ptranscend; 

theroots*NSolve[{r=-0, r==revfunc},{r,Q}]; 
theroots®N[{r,Q} /. theroots]; 

If[theroots=!={r,Q} && theroots !={},{ 
theroots=DeleteCases[theroots,{r_,Q_Complex /; 
Abs[Im[Q]]>.l}]; 

theroots=ReplaceAll[theroots,{r_,Q_} -> {r,Re[Q]}]; 
rmtranscend; 

Q=.; 

theroots=Table[FindRoot[revfunc*=0, 

{Q,theroots[[m,2]]}],{m,Length[theroots]}]; 
theroots*{0,Q} /. theroots; 
theroots*N[Union[theroots,{}]]; 

theroots=Union[theroots, theroots+ N[Table[{0,2Pi}, 

{i,Length[theroots]}]]]; 

theroots=Select[theroots,(#[[2]]<=N[upperlimit] && 

#[[2]] >= N[lowerlimit])&]; 

theroots=Union[N[(Round[100*theroots])/100],{}]; 

}].* 

rmtranscend; 

If[limits}cnown=!*y, { 

If[theroots=!={r,Q} && theroots!={} && 

Length[theroots]>=1,{ 

Print [" '• ] ; 
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Print[“The polar curve r is equal to zero when Q is: 

*•, Map [ topi, Transpose [ theroots ] [ [ 2 ] ] ] ] 

}].* 

Pause[8]; 

newcurve=Input["The polar curve is plotted from Q=o 

to Q-2 Pi. If you want to change the limits type y, otherwise 

type n."]; 

If[newcurve=®=y,{ 

lower1imit^Input["Input new lower 1imit.“]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi“>Pi; 

While[lowerlimit <-NtPi], {lowerlimit=Input[’'The lower limit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sgrt; 
lowerlimit=lowerlimit /. pi->Pi; 

}]; 

upperlimit=Input["Input new upper limit."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

}]; 

If [N[upperlimit]<N[lowerlimit] ,lowerlimit=N[lowerlimit-2Pi] ]; 

rayl=N[(func /. Q->lowerlimit)]; 
ray2=N[(func /. Q->upperlimit)]; 

llequivalent=N[Mod[lowerlimit, 2Pi]]; 
ulequivalent=N[Mod[upperlimit, 2Pi]]; 

If[lowerlimit < 0, llequivalent= 2 N[Pi]+lowerlimit, 
llequivalent= lowerlimit]; 

If[upperlimit < 0, ulequivalent= 2 N[Pi]+upperlimit, 
ulequivalent= upperlimit]; 

Q—.; 

If [ (rayl==ray2) && (llequivalent—ulequivalent) j j 
(rayl==0 && ray2==0),{ 

PolarPlot[func,(Q, topi[lower1imit], topi[upperlimit]}, 
AxesLabel -> {"x","y"}, AspectRatio -> Automatic, 

PlotRange-> {{-1.02 max, 1.02 max),{-1.02 max, 1.02 max)}]), 
{ 

Show[PolarPlot[func,{Q, topi[lowerlimit], topi[upperlimit]), 
DisplayFunction -> Identity, AxesLabel -> {"x","y"), 
AspectRatio -> Automatic, PlotRange-> {{-1.02 max,1.02 max), 
{-1.02 max, 1.02 max))], 

Graphics[{ 

Line[{{0,0),{Cos[lowerlimit]*rayl, Sin[lowerlimit]*rayl))], 


134 







Line[{ (0,0), {Cos[upperlimit]*ray2, Sin[upperliinit]*ray2))], 


)], 

DisplayFunction-> $DisplayFunction]}] 

)3 

)]; 

ptranscend; 

theroots=NSolve[{r==0, r==revfunc},{r,Q)]; 
theroots=N[(r,Q) /. theroots]; 

If[theroots=!={r,Q} && theroots !={},{ 
therootS’^DeleteCases [theroots, {r_,Q_Complex /; 
Abs[Iin[Q]]>.l}]; 

theroots=ReplaceAll[theroots,{r_,Q_) -> {r,Re[Q])]; 
mtranscend ; 

Q=. ; 

theroots=Table [ FindRoot [ revfunc“0, 

{Q, theroots [ [in/2] ]) ], {m, Length [theroots]} ]; 
theroots={0,Q) /. theroots; 
theroots*N[Union[theroots,()]]; 

theroots=Union[theroots, theroots+ N[Table[(0,2Pi), 

{i,Length[theroots])]]]; 

theroots=Select[theroots,(#[[2]]<=N[upperlimit] && 

#[[2]] >« N[lowerlimit])&]; 

theroots«Union[N[(Round[100*theroots])/100],{)]; 

}]? 

If[theroots ===(r,Q}||Length[theroots]<=1, { 
area^NIntegratei.5*(func)^2, {Q,topi[lowerlimit], 
topi[upperlimit])]; 
cltranscend; 

Print[" "]; 

Print["The limits of integration for the region bounded by 
the"]; 

Print["polar curve r = ",func, ", and the origin are:"]; 
Print[" "]; 

Print[" a=",topi[Chop[lowerlimit,.01]]," and b«", 

topi[Chop[upperlimit,.01]],"."]; 

Print[" "]; 

Print["The area bounded by the polar curve and the origin is 

It 

N[area,3]])]; 

If[theroots =!=(r,Q) && theroots!={) && Length[theroots]>1, ( 

theroots=Transpose[theroots]; 
theroots=theroots[[2]]; 
theroots=Sort[Union[ReplaceAll[Join[ 

{N[(Round[l00*lowerlimit])/100]), 

theroots,{N[(Round[100*upperlimit])/100]}],0.->0],{}]]; 
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rmtranscend; 

For[i=l, i<=*Ijength[theroots]-l, {ival=i; ++i; area=.; 
Q=(theroots[[ival]] + theroots[[ival+1]])/2; 

If[func=!=revfunc, { 

If[N[func] >0, { Q».; 

area*NIntegrate[.5*(func)*2,{Q,topi[theroots[[ival]] ], 
topi[theroots[[ival+1]] ])]; 

AppendTo[totalarea, Abs[area]]; 

AppendTo[regions, ival] }] 

}]; 

If[func===revfunc, { 

If[Abs[N[func]] >0, { Q=.; 

area=»NIntegrate[ .5* (func) ^2, {Q,topi[theroots[ [ival] ] ], 
topi[theroots[[ival+1]] ]}]; 

AppendTo[totalarea , Abs[area]]; 

AppendTo[regions , ival] }] 

)] 

}]; 

area=Apply[Plus, totalarea]; 
nuinregions=Length[regions]; 
cltranscend; 

If[numregions==l, {Q=.; 

Print [" '* ]; 

Print["The limits of integration for the region bounded by 
the"]; 

Print["polar curve r * ", func," and the origin are;"]; 

Print[" "]; 

Print[" a=",topi[Chop[theroots[[ regions[[l]] 

]],.01]], 

" and b=",topi[Chop[theroots[[regions[[1]] + 1]],.01]],"." ]; 
Print[" "]; 

Print["The area bounded by the polar curve and the origin 
is ",N[area,3] ])]; 

If[numregions >1, { 

Print[" "]; 

Print["The limits of integration for the regions which 
are bounded"]; Q=.; cltranscend; 

Print["by the polar curve r = ", func," and the 
origin are:"]; Print[" "]; 

If[Abs[N[Mod[theroots[[regions[[1]] ]]+2Pi,2Pi]] - 
N[Mod[theroots[[regions[[numregions]]+l]]+2Pi,2Pi]]]<.l,{ 

If[theroots[[regions[[1]] ]]<0 && 
theroots[[regions[[numregions]] ]]>=N[Pi], 

Print[" al=", 

topi[Chop[N[theroots[[regions[[numregions]] 
]]-2Pi],.01]],"bl=", 

topi[Chop[theroots[[regions[[1]]+l]],.01]] ]]; 
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If[therootsC[regions[[1]] ]]<0 ti 
theroots[[regions[[nunregions]] ]]<«N[Pi], 

Print[" al*", 

topi[Chop[theroots[[regions[[numregions]] ]],.01]], 

" and bl-'‘,topi[Chop[N[theroots[(regions[[l]]+l]]+2Pi],.01]] 

]]; 

If[theroots[[regions[[1]] ]]>*0 && 
theroots[[regions[[nunregions]] ]]>0, 

Print[" al®", 

topi[Chop[N[theroots[[regions[[nunregions]] ]]-2Pi],.01]], 

" and bl='*,topi[Chop[theroots[ [regions[ [1] ]+l] ], .01] ] ]]; 

For[j=3,j<*numregions, {jval=j; ++j; 

Print[" a",jval-1, 

topi[Chop[theroots[[regions[[jval-l]] ]],.01]], 

” and b”,jval-l,•*=•*, 

topi[Chop[theroots[[regions[[jval-l]]+!]],.01]] ])] 

), 

{For[j=l, j<=nunregions,{jval=j; ++j; 

Print[" a", jval,•*=•*, 

topi[Chop[theroots[[ regions[[jval]] ]],.01]], 

*' and b**,jval,•*='', 

topi[Chop[theroots[[ regions[[jval]]+l ]],.01]] ] )]; 

)]; 

Print[" "]; Q*.; 

Print["The total area bounded by the polar curve and the 
origin is"]; 

Print[" ",N[area,3]]; 

}]; 

)]; 

] 


ptranscend:=Module[(), 

sin[X_]:= x - X^3/6 + X^5/120 - X^7/5040 + X^9/362880 - 

X^ll/39916800+X^13/6227020800-X^15/1307674368000+ 
X^17/355687428096000-X^19/121645100408832000; 
X^21/51090942171709440000 - 
X^23/25852016738884976640000 + 
X^25/15511210043330985984000000; 
cos[x_]:= 1 - x^2/2 + X^4/24 - x^6/720 + x^8/40320 - 
X^lO/3628800 + X^12/479001600; 
tan[x_]:= x + X^3/3 + 2x^5/15 + 17x^7/315 + 62x^9/2835 + 
1382x^11/155925; 

COt[X_]:= 1/x -x/3 -X''3/45 - 2x''5/945 - X''7/4725 - 
2x^9/93555; 

csc[x_]:= 1/x + x/6 + 7x^3/360 + 31x^5/15120 + 

127x^7/604800 + 73x^9/3421440; 
sec[x_];= 1 + x^2/2 + 5x^4/24 + 61x^6/720 + 277x^8/8064 + 
50521x^10/3628800; 
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exp[x_]:* 1 + X +X^2/2 + X^3/6 + X''4/24 + X*5/120 + X''6/720 + 
X^7/5040 + X''8/40320 + X^9/362880] 

cltranscend:=Modulet{),sin[x_]=.; cos[x_]=.; tan[x_]=.; 
cot[x_3».; csc[x_]».; sec[x_]*.; exp[x_]».] 

rmtranscend:^Module[{), sin[x_]:=Sin[x]; cos[x_]:=Cos[x]; 
tan[x_]:=Tan[x]; cot[x_];=Cot[x]; csc[x_]:=Csc[x]; 
sec[x_]:=Sec[x]; exptx_]:=Exp[x] ] 

topic-.26]:*-Pi/12; topi[-.28]:=-Pi/ll; topi[-.29]:*-Pi/ll; 
topic-.31]:=-Pi/10; topic-.34]:=-Pi/9; topiC-.35]:=-Pi/9; 
topic-.39]:*-Pi/8; topiC-.44]:=-Pi/7; topiC-.45]:=-Pi/7; 
topic-.52]:=-Pi/6; topi{-.57]:=-2Pi/ll;topiC-.52];=-Pi/5; 
topi{-.63]:=-Pi/5; topi{-.69]:=-2Pi/9; topi{-.70]:=-2Pi/9; 
topiC-.78]:=-Pi/4; topiC-.79]:=-Pi/4; topiC-.85]:=-3Pi/llj 
topiC-.86] ;»-3Pi/ll; topic-.94] :=-3Pi/10; topiC-1.04]:=-Pi/3; 
topic-1.05]:=-Pi/3; topiC-1.14]:=-4Pi/ll; topiC-1.17]:=-3Pi/8; 
topi{-l.18]:=-3Pi/8; topiC-1.25]:=-2Pi/5; topi{-l.26];=-2Pi/5; 
topiC-1.30]:=-5Pi/12; topiC-1.31]:=-5Pi/12; 

topiC-1.34]:®-3Pi/7; topiC-1.35]:=-3Pi/7; topiC-1.39]:=-4Pi/9; 
topiC-1.40] :=-4Pi/9; topiC-1.42] :*-5Pi/ll; 

topiC-1.43]:=-5Pi/ll; topiC-1.57]:=-Pi/2; 

topiC-1.7l];=-6Pi/ll; topiC-1.74]:=-5Pi/9; 

topiC-1.75];=-5Pi/9;topiC-1.79]:=-4Pi/7;topiC-1.80]:=-4Pi/7; 
topiC-1.83];=-7Pi/12; 

topiC-1.88] :=-3Pi/5; topiC-1.96] :=-5Pi/8; 

topiC-1.99]:=-7Pi/ll; topiC-2.00]:=-7Pi/ll; 

topic-2.09] :=-2Pi/3; topic-2. 19] -7Pi/10; 

topiC-2.20]:=-7Pi/10; topiC-2.24]:=-5Pi/7; 

topiC-2.28]:=-8Pi/ll; topiC-2.35]:=-3Pi/4; 

topiC-2.36]:*-3Pi/4; topiC-2.44]i=-7Pi/9 ? topiC-2.51]:=-4Pi/5; 
topiC-2.57] :=-9Pi/ll; topiC-2.61] ;=-5Pi/6; 

topiC-2.62]:=-5Pi/6; topiC-2.69]:=-6Pi/7; topiC-2.74]:=-7Pi/8; 
topiC-2.75]:=-7Pi/8; topiC-2.79]:=-8Pi/9; 

topiC-2.82]:=-9Pi/10; topiC“2.83]:=-9Pi/10; 

topic-2.85] :=-10 Pi/ll; topic-2.86] :=-lOPi/ll; 
topiC-2.87]:=-llPi/12; topic-2.88]:=-llPi/12; 
topiC-3.14]:=-Pi; topic.26]:=Pi/12; topiC.28];=Pi/ll; 
topi{.29]:=Pi/ll7 topic.31]:=Pi/10; topiC.34]:=Pi/9; 
topic.35];=Pi/9; topiC.39]:=Pi/8; topiC.44]:=Pi/7; 
topic *45] :=Pi/7;topic .52] :=Pi/6;topic *57] :=2Pi/ll;topic . 62 ] : = 
Pi/5; topic.63]:=Pi/5; topic.69];*2Pi/9; topic.70]:=2Pi/9; 
topic.78];=Pi/4; topic.79]:=Pi/4; topiC.85]:=3Pi/ll; 
topic.86]:=3Pi/ll; topic.94];=3Pi/10; topiC1.04];=Pi/3; 
topi{l.05]:=Pi/3; topiCl.14]:=4Pi/ll; topifl.l7];=3Pi/8; 
topiC1.18]:=3Pi/8; topiCl.25]:=2Pi/5; topiCl.26]:=2Pi/5; 
topiC1.30];=5Pi/12; topiC1.31]:=5Pi/12; topiC1.34]:=3Pi/7; 
topifl.35]:=3Pi/7; topiC1.39];=4Pi/9; topiC1.40]:=4Pi/9; 
topiCl.42]:=5Pi/ll; topiCl.43];=5Pi/ll; topicl.57]:=Pi/2; 
topiC1.71]:=6Pi/ll; topiCl.74];=5Pi/9; topiC1.75];=5Pi/9; 
topiCl.79];=4Pi/7; topiCl.80]:=4Pi/7; topiCl.83]:=7Pi/12; 
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topi[1.88]:»3Pi/5; topi[l.96]:-5Pi/8; topi[1.99]:=7Pi/ll; 
topi[2.00]:»7Pi/ll; topi[2.09]:»2Pi/3; topi[2.19];« 7Pi/lO; 
topi[2.20]:*7Pi/10; topi[2.24]:«5Pi/7; topi[2.28]:«8Pi/ll; 
topi[2.35];*3Pi/4; topi[2.36]i-3Pi/4; topi[2.44]:*7Pi/9; 
topi[2.51]:»4Pi/5; topi[2.57]:»9Pi/ll; topit2.61]:=5Pi/6; 
topit2.62];*5Pi/6; topi[2.69];=6Pi/7; topi[2.74]:=7Pi/8; 
topi[2.75]:=7Pi/8; topi[2.79];=8Pi/9; topif2.82]:*9Pi/10; 
topi[2.83]:=9Pi/10; topi[2.85]:-10Pi/ll;topi[2.86]:=10Pi/ll; 
topi[2.87];*llPi/l2; topi[2.88]:sllPi/12; topi[3.14]:*Pi; 
topi[3.40]:»13Pi/12;topi[3.42]:*12Pi/ll;topi[3.43]:-12Pi/ll; 
topi [3.45]:»llPi/10; topi [3.46] :«llPi/10;topi[3.49] :-10Pi/9 ; 
topi[3.59]:-8Pi/7; topi[3.53]:=9Pi/8; topi[3.66]:«7Pi/6; 
topi[3.67] :=7Pi/6,* topi[3.71] :-13Pi/ll; topi[3.76]:-6Pi/5; 
topi[3.77]:=6Pi/5; topi[3.83]:«llPi/9; topi[3.84]:»llPi/9; 
topi[3.92]:*5Pi/4; topi[3.93]:*5Pi/4; topi[3.99]:»14Pi/ll; 
topi[4.00]:=14Pi/ll; topi[4.03]:=9Pi/7; topi[4.04];*9Pi/7; 
topi[4.08]:=13Pi/10; topi[4.18]:=4Pi/3; topi[4.19]:=4Pi/3; 
topi[4.28]:=15Pi/ll; topi[4.31]:=llPi/8; topi[4.32]:=llPi/8; 
topi[4.39]:=7Pi/5; topi[4.40]:=7Pi/5; topi[4.45]:=17Pi/l2; 
topi[4.48];=10Pi/7; topi[4.49]:=10Pi/7; topi[4.53]:=13Pi/9; 
topi[4.54]:=13Pi/9; topi[4.56]:=16Pi/ll;topi[4.57]:*16Pi/ll; 
topi[4.71]:=3Pi/2; topi[4.85]:=17Pi/ll; topi[4.86]:*17Pi/ll; 
topi[4.88]:®14Pi/9; topi[4.89]:=14Pi/9; topi[4.93]:=llPi/7; 
topi[4.94]:=llPi/7; topi[4.97]:=19Pi/12; topi[5.02]:=8Pi/5; 
topi[5.03]:=8Pi/5; topi[5.10]:=13Pi/8; topi[5.11]:=13Pi/8; 
topi[5.14]:=18Pi/ll; topi[5.23]:=5Pi/3; topi[5.24]:=5Pi/3; 
topi[5.34]:=17Pi/10; topi[5.38]:*12Pi/7; topi[5.39]:=12Pi/7; 
topi[5.42]:*19Pi/ll; topi[5.43]:=19Pi/ll; topi[5.49]:=7Pi/4; 
topi[5.50]:=7Pi/4; topi[5.58]:=16Pi/9; topi[5.59]:=16Pi/9; 
topi[5.65]:=9Pi/5; topi[5.71];=20Pi/ll; topi[5.75]:=llPi/6; 
topi[5.76]:=llPi/6; topi[5.83]:=13Pi/7: topi[5.89]:=15Pi/8; 
topi[5.93]:=17Pi/9; topi[5.96]:=19Pi/10;topi[5.97]:=19Pi/10; 
topi[5.99]:=21Pi/ll;topi[6.00]:=21Pi/ll;topi[6.02]:=23Pi/12; 
topi[6.28]:=2Pi; topi[6.54]:=25Pi/12; topi[6.56]:*23Pi/ll; 
topi[6.57]:=23Pi/ll;topi[6.59]:=21Pi/10;topi[6.60]:=21Pi/10; 
topi[6.63]:=19Pi/9; topi[6.67];*17Pi/8; topi[6.68]:=17Pi/8; 
topi[6.73]:=15Pi/7; topi[6.80]:=13Pi/6; topi[6.81]:=13Pi/6; 
topi[6.85]:=24Pi/ll; topi[6.91]:=llPi/5; topi[6.98]:=20Pi/9; 
topi[7.06]:=9Pi/4; topi[7.07];=9Pi/4; topi[7.l3]:=25Pi/ll; 
topi[7.14]:=25Pi/ll; topi[7.18]:=16Pi/7;topi[7.22]:=23Pi/10; 
topi[7.23]:=23Pi/10; topi[7.33]:=7Pi/3; topi[7.42]:=26Pi/ll; 
topi[7.43]:=26Pi/ll; topi[7.46];=19Pi/8; topi[7.53]:=12Pi/5; 
topi[7.54]:=12Pi/5; topi[7.59]:=29Pi/12; topi[7.62]:=17Pi/7; 
topi[7.63]:=17Pi/7; topi[7.67]:=22Pi/9; topi[7.68]:=22Pi/9; 
topi[7.7l]:=27Pi/ll; topi[7.85]:=5Pi/2; topi[7.99]:=28Pi/ll; 
topi[8.00]:=28Pi/ll; topi[8.02];=23Pi/9; topi[8.03]:=23Pi/9; 
topi[8.07]:=l8Pi/7; topi[8.08]:=18Pi/7; topi[8.11]:=31Pi/12; 
topi[8.12]:=31Pi/12; topi[8.16]:=13Pi/5; topi[8.17]:=13Pi/5; 
topi[8.24]:=21Pi/8; topi[8.25]:=21Pi/8; topi[8.28];=29Pi/ll; 
topi[8.37]:=8Pi/3; topi[8.38]:=8Pi/3; topi[8.48]:=27Pi/lO; 
topi[8.52]:=i9Pi/7; topi[8.53]:=19Pi/7; topi[8.56]:=30Pi/ll; 
topi[8.57]:=30Pi/ll; topi[8.63]:*llPi/4; topi[8.64]:=llPi/4; 
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topi[8.72] 
topi[8.80] 
topi[8.97] 
topi[9.07] 
topi[9.13] 
topi[9.42] 


-25Pi/9; topi[8.73]:=25Pi/9; topi[8.79]:*14Pi/5; 
»14Pi/5; topi[8.85]:=31Pi/ll; topi[8.90];=17Pi/6; 
=20Pi/7; topit8.98]:=20Pi/7; topi[9.03];»23Pi/8; 
=26Pi/9; topi[9.08]:=26Pi/9; topi[9.11];*29Pi/10; 
=32Pi/ll;topi[9.14]:=32Pi/ll;topi[9.16]:*35Pi/12; 
=3Pi; topi[x_]:=x; 
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yix. VILB POIARl.MA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 
polar2:-Module[{ 

), 

cltranscend; 

fund-input["Enter one of the polar curves which bounds the 
region. 

r = " ]; 

fund-fund /. sqrt->Sqrt; 
fund-fund /. pi->Pi; 
fund-fund /. q->Q; 

func2=Input["Enter the other polar curve which bounds the 
region. 

r = " ]; 

func2=func2 /. sqrt->Sqrt; 
func2=func2 /. pi“>Pi; 
func2=func2 /. q->Q; 

liinitsknown=Input["If the limits are known type y, 
otherwise type n."]; 

If[limitsknown—y, { 

lowerlimit-Input["Input lower limit."]; 
lowerlimit-lowerlimit /. sqrt->Sqrt; 
lowerlimit-lowerlimit /. pi->Pi; 
upper1imit-Input["Input upper 1imit."]; 
upperlimit-upperlimit /. sqrt->Sqrt; 
upperlimit-upperlimit /. pi->Pi; 

If[N[upperlimit]<N[lowerlimit],lowerlimit=N[lowerlimit-2Pi]] 

), 

{lowerlimit-0; upperlimit=2 N[Pi])]; 
bothfunc- {fund, func2} ; 
ptranscend; 

If [{fund, func2 ) = ! = ({ fund, func2) /. Sqrt[x_]->x), 

theroots=NRoots[funcl^2—func2''2 ,Q], 
theroots=NRoots[funcl—func2, Q]; 

]? 


If[theroots-!-False, { 

theroots- Table[theroots[[i,2]],{i,Length[theroots]}]; 
theroots=DeleteCases[theroots, Q_Complex /; Abs[Im[Q]]>.1]; 
theroots-ReplaceAll[theroots, Q_ -> Re[Q] ]; 
rmtranscend; 

Q-. ; 
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theroots«Table[FindRoot[funcl*»func2# {Q,theroots[[m]]) ], 

{m, Length[theroots])]; 

theroots- Q /. theroots; 

theroots»N[Union[theroots , {>31? 

theroots«Union[theroots, theroots + N[Table[2Pi, 

{j, Length[theroots]}]]]; 

testroots^Union[Select[N[Round[theroots*100]/100], 
(#>=N[-Pi] && #<-N[(5/2) Pi])4],()]; 
theroots^Select[theroots,(#<» N[upperlimit] && 

#>= N[lowerlimit])i]; 

theroots»Union[N[(Round[100*theroots]l/lOO],{)]; 
theroots-Map[topi, theroots]; 
testrootSBMap[topi, testroots]; 

)]; 

theroots^Transpose[Select[Tablet(theroots[[i]3 > 
(N[Abs[funri-func23 /. Q-> theroots[[i]] ])<.05), 

{i,1,Length[theroots])],(#[[2]3===True)&]][[1]]; 

If[theroots==={)[[1]],theroots=()]; 

max={0,0}; 

For[rr=l, rr<=2, {rval=rr; ++rr; 

max[[rval]]-Max[OeleteCases[Table[bothfunc[[rval]] 

/. Q->values, (values, 0, 2Pi, .25}], x_Complex]] )]; 
max=Max[max]; 

plotroutine; 

ptranscend; 

For[hh=l, hh<=2, (hval=hh; ++hh; 

otherroots=NSolve[(r==0, r==bothfunc[[hval]]), {r,Q)]; 
otherroots=N[(r,Q) /. otherroots]; 


If[otherroots=! = {r,Q} & & otherroots ! = {},( 
otherroots=DeleteCases[otherroots,(r_,Q_Complex /; 
Abs[Im[Q]]>.l)]; 

otherroots^ReplaceAl1[otherroots,{r_,Q_) -> (r,Re[Q]}]; 
rmtranscend; 

Q=. ; 


otherroots=Table[FindRoot[bothfunc[[hval]]==0, 

(Q,otherroots[[m,2]]}],(m,Length[otherroots])]; 
otherroots=(0,Re[Q]) /. otherroots; 
otherroots=N[Union[otherroots,()]]; 

otherroots=Union[otherroots, otherroots+ N[Table[{0,2Pi), 
(i,Length[otherroots]}]]]; 

otherroots=Select[otherroots,{#[[2]]<=N[upperlimit] && 
#[[2]] >= N[lowerlimit])&]; 

otherroots=Union[N[(Round[100*otherroots])/100],()]; 

}]? 
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otherroots=UniontReplaceAll[otherroots,{0,6.28}->{0,0)l,{)]; 
otherroots»Union[ReplaceAll[otherroots, (Q,r)->{{Q,r))],{)]; 
If[hval*=l, funclroots^otharroots, func2roots®otherroots]; 
}]; 


allroots={funclroots, func2roots}; 
endprogram^O; 
runoption*0; 
diffcheck; 

If[(theroots===False [j testdiff===yes) && runoption==0 && 
Length[funclroots]<=l && Length[func2roots]<=1,{ 

If[liinitsknown===n, { 

If[Length[theroots]==1/ 

Print[" 

”The two polar curves intersect at Q = ", 
topi[theroots[[1]] ] ]]; 

If[Length[theroots] >1, 

Print["The two polar curves intersect at Q: "," 

", Map[topi, theroots] ]]; 

Pause[8]; 

newcurve=Input["The polar curves are plotted from Q=0 
to Q=2Pi. If you want to change the limits type y, otherwise 
type n."]; 

If[newcurve===y,{ 

lowerlimit=Input["Input new lower limit."]; 
lowarlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

While[lowerlimit <-N[Pi], {lowerlimit=Input["The lower limit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

)]? 

upperlimit=Input["Input new upper limit."]; 
upperlimit-upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit*upperlimit /. pi->Pi; 

)]? 

If [N[upperlimit ]<N[lowerl i t] ,lowerlimit=N[lowerlimit-2Pi] ]; 
rmtranscend; plotroutine, 

)V 

}]; 

complexcheck; negcheck; 

If[Length[complexfuncl]<=2 && Length[complexfunc2]<=2 && 
endprogram==0, findarea2]; 

If[(Length[complexfunci] >2 j| Length[complexfunc2] >2) && 
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endprograin-=0, findareal]; runoption=l 

)]; 


If[testdiff==*yes && (Length[fvnclroots] >1 || 

Length[func2roots]>l) && runoption==0, { 

If[limitsknown==*n, { 

If[Length[theroots]==1, 

Print [•• 

"The two polar curves intersect at Q = 
topi[theroots[[1]] ] ]]; 

If[Length[theroots] >1, 

Print["The two polar curves intersect at Q: "," 

", Map[topi, theroots] ]]; 

Pause[8]; 

newcurve=Input["The polar curves are plotted from Q=0 
to Q=2Pi. If you want to change the limits type y, otherwise 
type n."]; 

If[newcurve===y,{ 

lowerlimit=Input["Input new lower limit."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

While[lowerlimit <-N[Pi], {lowerlimit=Input["The lower limit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

}]? 

upperlimit=Input["Input new upper limit."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

)]; 

If [N[upperlimit]<N[lowerlimit],lowerlimit=N[lowerlimit-2Pi] ]; 
rmtranscend; plotroutine; 

)]; 

)]; 

complexcheck; negcheck; 

If[Length[complexfuncl]<=2 && Length[complexfunc2]<=2 && 
endprogram==0, findarea2]; 

If[(Length[complexfuncl] >2 || Length[complexfunc2] >2) && 
endprogram==0, findareal]; runoption=l 
}]? 

If[(theroots===False j[ Length[theroots] <®1) && 
runoption==0 && 

( Length[funclroots] >1 |i Length[func2roots]>1), { 

If[limitsknown===n, { 
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If[Length[theroots]**1, 

Print [" •*," 

The two polar curves intersect at Q * ", 
topi[theroots[[1]] ] ]]; 

Pause[8]; 

newcuzn^e-InputCThe polar curves are plotted from Q=0 
to Q>:2Pi. If you want to change the limits type y, otherwise 
type n. •*]; 

If[newcurve===y,{ 

lowerlimit=Input["Input new lower limit."]; 
lowerlimit»lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

While[lowerlimit <-N[Pi], {lowerlimit=Input["The lower limit 
can not be less than -Pi» input lower limit again."]; 
lowerlimit-lowerlimit /. sqrt->Sqrt; 
lowerlimit-lowerlimit /. pi->Pi; 

)]; 

upperlimit=Input["Input new upper limit."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi“>Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

)]; 

If [N [upperlimit ]<N[ lower limit] ,lowerlimit=N[lowerlimit-2Pi] ]; 
rmtranscend; plotroutine; 

)]; 

)]; 

complexcheck; negcheck; 

If[Length[complexfunci]<=2 && Length[complexfunc2]<=2 && 
endprogram==0 , findarea2]; 

If[(Length[complexfunci] >2 |j Length[complexfunc2] >2) && 
endprogram==0, findareal];runoption*l 
}]; 

control=0; 

If[Length[theroots]>=2 && limitsknown===n && runoption==0,{ 
Print["The two polar curves intersect at Q: "," 

", Map[topi, theroots] ]; 

Pause[8]; 

newcurve=Input["The polar curves are plotted from Q=0 
to Q=2Pi. If you want to change the limits type y, otherwise 
type n."]; 

If[newcurve===y,( 

lowerlimit=Input["Input new lower limit."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

While[lowerlimit <-N[Pi], (lowerlimit=Input["The lower limit 
can not be less than -Pi, input lower limit again."]; 
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loverllmit«lowerlimit /. 8qrt->Sqrt; 
lowerlinltslowerliinit /. pi->Pi; 

)]; 

upperlliDit->Input[''Input new upper limit."]; 
upperlimit-upperlimit /. sqrt->Sqrt; 
upperlimitsupperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit>Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit-upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

)]; 

If [N[ upperl imit] <N[ lower limit] ,lowerlimitsN[lowerlimit-2Pi] ]; 

)]; 

badlimitssSelect[N[theroots], 

(#>N[lowerlimit] && #<N[upperlimit])&]; 

While[Length[badlimits]>0, { 

If[control==l, { 

Print["For Q=-Pi to Q=5/2 Pi the two polar curves intersect 
at Q:","", Map[topi, testroots] ]; Pause[8]}]; 
control®control+l; 

lowerlimit=Input["A point of intersection can not lie 

between the two limits. Input new lower limit."]; 

lowerlimit=lowerlimit /. sqrt->Sqrt; 

lowerlimit=lowerlimit /. pi~>Pi; 

upperlimit=Input["Input new upper limit."]; 

upperlimit=upperlimit /. sqrt->Sqrt; 

upperlimit=upperlimit /. pi->Pi; 

If [N[upperlimit]<N[lowerlimit] ,lowerlimit=N[lowerlimit-2Pi] ]; 
badlimits=Select[N[testroots], 

(#>N[lowerlimit] && #<N[upperlimit])&]; 

}]; 

rmtranscend; plotroutine; complexcheck; negcheck; 

If[(Length[complexfuncl]<=2 && Length[complexfunc2]<=2) && 
endprogram==0,findarea2]; 

If[(Length[complexfunci] >2 |j Length[complexfunc2] >2) && 
endprogram==0,findareal]; runoption=l 
)]; 

onceprint=0; 

If[Length[theroots]>=l && 1imitsknown===y && runoption==0,( 
badlimits=Select[N[theroots], 

(#>N[lowerlimit] && #<N[upperlimit])&]; 

While[Length[badlimits]>0, ( 

If[onceprint==0,( 

Print["The two polar curves intersect at Q: "," 

", Map[topi,theroots] ]; onceprint*l)]; 

Pause[8]; 

If[control==l, I 

Print["For Q=-Pi to Q=5/2 Pi the two polar curves intersect 
at Q:","", Map[topi, testroots] ]; Pause[8])]; 
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control=control+l; 

lowerlimlt»Input['*A point of intersection can not lie 
between the two limits. Input new lower limit.**]; 
lowerlimit^^lowerlimit /. sqrt->Sqrt; 
lowerlimit^lowerlimit /. pi->Pi; 
upperlimit=Input[**Input new upper limit.**]; 
upperlimit»upperlimit /. sqrt->Sqrt; 
upperlimit^upperlimit /. pi->Pi; 

If [N[upperlimit]<N[lowerlimit] ,lowerlimit=N[lowerlimit-2Pi] ]; 
badlimitssSelect[N[testroots], 

(#>N[lowerlimit] 6& #<N[upperlimit])4]; 

If[Length[badlimits]==0/ { 
rmtranscend; plotroutine}]; 

}]; 

complexcheck; negcheck; 

If[(Length[complexfuncl]<=2 && Length[complexfunc2]<=2) && 
endprogram==0,findarea2]; 

If[(Length[complexfuncl] >2 {| Length[complexfunc2] >2) && 
endprogram==0 , findareal]; runoption=l 

)]; 

]; 


complexcheck:=Module [(}, 

complexfuncl=DeleteCases[DeleteCases[Table[N[fund /. 
Q->values], (values, N[lower1imit], N[upper1imit], .1}], 
r_Integer] , r_Real]; 

complexfunc2»=0eleteCases [DeleteCases [Table [N [ func2 /. 
Q->values], (values,N[lowerlimit],N[upperlimit], .1)], 
r_Integer], r__Real] ]; 

diffcheck:=Module[{}, 

difference=DeleteCases[Table[funcl-func2 /.Q->values, 
(values, lowerlimit+.Ol, upperlimit-.Ol,.25}], x_Complex]; 
If[Max[difference]<=0 && Min[difference]<0 {j 
Max[difference]>0 && Min[difference]>=0, 
testdiff=yes, testdiff=no] ]; 

findareal:=Module[{), 

For[tt=l, tt<=2, {tval=tt; ++tt; 

If[allroots[[tval]] =1= {{Q,r)), { 
otherroots=Transpose[ allroots[[tva1]] ]; 
otherroots=otherroots[[2]]; 

otherroots=Sort[Union[ReplaceAll[Join[{N[lowerlimit]}, 
otherroots, {N[upperlimit]) ], 0.->0],n]]»’ 
otherroots=Select[otherroots,(#>=N[lowerlimit] && 
#<=N[upperlimit])&]; 

),{otherroots=ReplaceAll[{N[lowerlimit],N[upperlimit]}, 
0.->0]}]; 

rmtranscend; area=.; totalarea={); 

For[vv=l, vv<=Length[otherroots]-l, {vval=vv; ++vv; 

Q= (otherroots[ [wal] ] + otherroots[ [wal +l]])/2; 

If[ N[bothfunc[[tval]] ]>0, {Q=.; 
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area* NIntegrate[.5* bothfunct[tval]]^2, {Q, 

topi [otherroots [ [ wal ] ] ], topi [otherroots [ [ wal+1] ] ]} ]; 

AppendTo[totalarea, Abs[area] ]}]; 

)]; 

area» Apply[Plus, totalarea]; 

If[tval==l, areal=area]; 

If[tval»*2, area2»area]; 

)]; 

areas area2 - areal; cltranscend;Q*.; 
printfuncl ]; 

findarea2:^Module[{), 

Q=.; rmtranscend; area=NIntegrate[ .5*(func2''2 - funcl^2), 
{Q, topi[loverliinit], topi[upperliinit]}]; printfunc2]; 

printfuncl:=Module[{}, 

lf[area<=0, {tempsfuncl; funcl=func2; func2steiDp} ]; 

Print [•» ••]; 

Print[" "]; 

Print['*The outer polar curve is r = ", func2]; 

Print ["The inner polar curve is r * ", fund]; 

Print["The area of the region bounded by the two"]; 

Print["polar curves is ",Abs[area]] ]; 

printfunc2;=Module[{), 
cltranscend; 

If[area<sO, {tempsfuncl; funcl==func2; func2=temp) ]; 

Print[" "]; 

Print[" "]; 

Print["The outer polar curve is r * ", func2]; 

Print ["The inner polar curve is r * ", fund]; 

Print["The limits of integration are a=", 
topi[N[Round[lowerlimit*100]/l00]], 

" and b=", topi[N[Round[upperlimit*100]/100]] ]; 

Print["The area of the region bounded by the two"]; 

Print["polar curves is ", Abs[area]] ]; 

negcheck:s Module[{), 

negfuncl=Min[DeleteCases[Table[fund /. Q->values, 

{values, lowerlimit+.01, upperlimit-.01, .25)], x_Complex]] 
negfunc2=Min[DeleteCases[Table[func2 /. Q->values, 

{values, lowerlimit+.01, upperlimit-.01, .25)], x_Complex]] 
lf[negfuncl<0 i| negfunc2<0,{cltranscend; 
lf[negfuncl<0 && negfunc2<0, 

{Print["Both functions have negative values for r."]; 
Print["This program can not solve such problems."]), 

If [negfuncKO, 

{Print["The function r = ",funcl," has negative values."] 
Print["This program can not solve such problems."] ), 
{Print["The function r = ",func2," has negative values."] 
Print["This program can not solve such problems."])] 

]; 
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endprogram^l 

}] 

]; 


plotroutine:* 

Module[{}, 

plotchunk=PolarPlot[ {fund, func2), {Q, N[loverliiDit], 

N[upper1imit]}, DisplayFunction->Identity, AxesLabel-> 

{"x","y"}, AspectRatio->Automatic, PlotRange-> 

{{-1.02 max, 1.02 max),{-1.02 max, 1.02 max))]; 

rayl=Re[N[ (fund /. Q-> lowerlimit) ] ]; 
ray2=Re[N[ (fund /. Q-> upperlimit) ] ]; 
ray3=Re[N[(func2 /. Q-> lowerlimit)]]; 
ray4*Re[N[(func2 /. Q-> upperlimit)]]; 
llequivalent=N[Mod[lowerlimit, 2Pi]]; 
ulequivalent=N[Mod[upperlimit, 2Pi]]; 

If[lowerlimit < 0, llequivalent= 2 N[Pi]+lowerlimit, 
llequivalent=lowerlimit]; 

If[upperlimit < 0, ulequivalent= 2 N[Pi]+upperlimit, 
ulequivalent=upperlimit]; 

Q=.; 

If[(rayl==ray2 && llequivalent==ulequivalent && 
ray3“ray4) I { (rayl*=0 && ray2~0 && ray3=*0 && ray4~0), 
{Show[plotchunk, DisplayFunction-> $DisplayFunction]), 
{Show[plotchunk, Graphics[{ 

Line[{{Cos[lowerlimit]*rayl. Sin[lowerlimit]*rayl), 

{Cos[lowerlimit]*ray3, Sin[lowerlimit]*ray3))], 

Line[{{Cos[upperlimit]*ray2, Sin[upperlimit]*ray2), 

{Cos[upperlimit]*ray4. Sin[upperlimit]*ray4))], 

)], DisplayFunction-> $DisplayFunction] 

}]] 

ptranscend:=Module[{), 

Sin[x_];= x - X^3/6 + x^5/120 - x^7/5040 + X^9/362880 - 

X^ll/39916800+x^l3/6227020800-X^15/1307674368000+ 
X^17/355687428096000-x^l9/121645100408832000; 
cos[x_];= 1 - X^2/2 + X^4/24 - x^6/720 + x''8/40320 - 
X^lO/3628800 + X^12/479001600; 
tan[x_];= x + X''3/3 + 2x^5/15 + 17x^7/315 + 62x''9/2835 + 
1382x^11/155925; 

cot[x_]:= 1/x -x/3 -X''3/45 - 2x^^57945 - x^7/4725 - 
2x^9/93555; 

CSC[x_];= 1/x + X/6 + 7x^3/360 + 31x''5/15120 + 

127x^7/604800 + 73x^9/3421440; 
sec[X_];= 1 + x^2/2 + 5x^4/24 + 61x^6/720 + 277x^8/8064 + 
50521x^10/3628800; 
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exp[x_]:» 1 + X +X*2/2 + X'^a/e + X^4/24 + x^5/120 + 
X^7/5040 + X*8/40320 + X^9/362880] 


I 


X^6/720 + 
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VIII. FILE POLPLOTl.MA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 

polplotl:* 

Module[{ 

func, limitsknown, theroots, lowerlimit, 
upperlimit, newcurve, curvesknown, temp 

), 

Q~ • » 

cltranscend; 

funcs^Input["Enter the polar curve which bounds the region. 


r = "]; 

func=func /. sqrt->Sqrt; 
func=func /. pi->Pi; 
func=func /. q->Q; 

revfunc=Replace[func, Sqrt[x_] -> x]; 
limitskno%m=Input["If the limits are known type y, 
otherwise type n."]; 

If[1imitsknown===y, 

{lowerlimit=Input["Input lower limit."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi“>Pi; 

While[lowerlimit <-N[Pi], {lowerlimit=Input["The lowerlimit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

)]; 

upperlimit=Input["Input upper limit."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=lnput["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperliiT\it=upperliniit /. pi->Pi; 

)] 

), 

{lower1imit=0.; upperlimit® 2Pi}]; 
rmtranscend; 

max=Max[DeleteCases[Table[func /. Q -> values, 

{values, 0, 2Pi, .25)],x_Complex]]; 

rayl=N[(func /. Q->lowerlimit)]; 
ray2=N[(func /. Q->upperlimit)]; 




llequlvalent»>N[Mod[lowerlimit, 2Pi] ]; 
ulequivalent»N[Mod[upperllBit, 2Pi]]; 

If[lowerlimit < 0, llequivalent» 2 N[Pi]-t-lowerlimit, 
llequivalent» lowerlimit]; 

If[upperlimit < 0, ulequivalent* 2 N[Pi]+upperlimit, 
ulequivalent» upperlimit]; 

Q® •» 

If[(rayl“ray2) && {llequivalent==ulequivalent)|j 
(rayl==0 && ray2=»0),{ 

PolarPlot[func,{Q, lowerlimit, upperlimit), 

AxesLabel -> {"x","y'*), AspectRatio -> Automatic, 

PlotRange-> {{-1.02 max, 1.02 max),{-1.02 max, 1.02 max}}]}, 
{ 

Show[PolarPlot[func,{Q, lowerlimit, upperlimit}, 
DisplayFunction -> Identity, AxesLabel -> {'*x'*,"y"}, 
AspectRatio -> Automatic, PlotRange-> {{-1.02 max,1.02 max}, 
{-1.02 max, 1.02 max}}], 

Graphics[{ 

Line[{{0,0},{Cos[lowerlimit]*rayl, sin[lowerlimit]*rayl}}], 
Line[{{0,0},{Cos[upperlimit]*ray2. Sin[upperlimit]*ray2}}], 

)], 

DisplayFunction-> $DisplayFunction ]}]: 
ptranscend; 

theroots=NSolve[{r®*0, r==revfunc},{r,Q}]; 
theroots=N[{r,Q) /. theroots]; 

If[theroots=l={r,Q) && theroots !={},{ 
theroots=OeleteCases[theroots,{r_,Q_Complex /; 
Abs[Im[Q]]>.l}]; 

theroots=ReplaceAll[theroots,{r_,Q_} -> {r,Re[Q]}]; 
rmtranscend; 

Q=. J 

theroots=Table[FindRoot[revfunc==0, 

{Q,theroots[[m,2]]}],{m,Length[theroots ]}]; 
theroots={0,Q} /. theroots; 
theroots=N[Union[theroots,{}]]; 

theroots=Union[theroots, theroots+ N[Table[{0,2Pi}, 

{i,Length[theroots]}]]]; 

theroots=Select[theroots,(#[[2]]<=N[upperlimit] && 

#[[2]] >=N[lowerlimit])&]; 

theroots=Union[N[(Round[100*theroots])/100],{}]; 

)]; 

rmtranscend; 

If[limitslaiown=!=y, { 

If[theroots=!={r,Q} && theroots!={} 

Length[theroots]>*1,{ 

Print [" '•]; 
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Print["The polar curve r is equal to zero when Q is; 

"jMapCtopijTransposeCtheroots][[2]] ] ] 

)]; 

Pause[4]; 

newcurve=Input["The polar curve is plotted from Q=0 

to Q=2 Pi. If you want to change the limits type y, otherwise 

type n."]; 

If[newcurve===y,{ 

lower1imit*Input["Input new lower 1imit."]; 
lowerlimit-lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

While[lowerlimit <-N[Pi], {lowerlimit^Inputf’The lower limit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

}]? 

upperlimit=Input["Input new upper limit."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

}]; 

rayl=N[(func /. Q->lowerlimit)]; 
ray2=N[(func /. Q->upperlimit)]; 

llequivalent=N[Mod[lowerlimit, 2Pi]]; 
ulequivalent=N[Mod[upperlimit, 2Pi]]; 

If[lowerlimit < 0, llequivalent= 2 N[Pi]+lowerlimit, 
llequivalent= lowerlimit]; 

If[upperlimit < 0, ulequivalent= 2 N[Pi]+upperlimit, 
ulequivalent= upperlimit]; 

Q=. ; 

If[(rayl==ray2) && (llequivalent==ulequivalent)|j 
{rayl==0 && ray2“0), { 

PolarPlot[func,{Q, lowerlimit, upperlimit), 

AxesLabel -> {"x","y"}, AspectRatio -> Automatic, 

PlotRange-> {{-1.02 max, 1.02 max),{-1.02 max, 1.02 max})]}, 
{ 

Show[PolarPlot[func,{Q, lowerlimit, upperlimit), 
DisplayFunction -> Identity, AxesLabel -> {"x","y"), 
AspectRatio -> Automatic, PlotRange-> {{-1.02 max,1.02 max), 
{-1.02 max, 1.02 max})], 

Graphics[{ 

Line[{{0,0),{Cos[lowerlimit]*rayl. Sin[lowerlimit]*rayl))], 
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Line[{{0,0),{Cos[upperlimit)*ray2, Sin[upperlimit]*ray2))), 


) 3 , 


DisplayFunction-> $DisplayFunction])] 
)]? 

)]; 

] 
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IX. riLB POLPLOT2.1IA 


(* Copyright Dennis A. Polaski, May 4,1993 *) 

polplot2:“Module[{ ), 

cltranscend; 

funcl=Input["Enter one of the polar curves which bounds the 
region. 

r = " ]; 

funcl=funcl /. sqrt->Sqrt; 
funcl=funcl /. pi->Pi; 
funcl=funcl /. q->Q; 

func2“Input["Enter the other polar curve which bounds the 
region. 

r = " ] ; 

func2=func2 /. sqrt->Sqrt; 
func2=func2 /. pi->Pi; 
func2=func2 /. q->Q; 

liinitsknown=Input[ "If the limits are known type y, 
otherwise type n."]; 

If[limitsknown===y, { 

lowerlimit“Input["Input lower limit."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit*lowerlimit /. pi->Pi; 
upper1imit“Input["Input upper 1imit."]; 
upperlimit“upperlimit /. sqrt->Sqrt; 
upperlimit“upperlimit /. pi->Pi), 

{lowerlimit=0; upperlimit=2 N[Pi]}]; 
bothfunc= {fund, func2); 
ptranscend; 

If [{fund, func2 } = ! = (( fund, func2} /. Sqrt[x_]->x), 

theroots=NRoots[func1^ 2==func2 ^ 2,Q], 
theroots=NRoots[funcl==func2, Q]; 

]; 


If[theroots=!“False, { 

theroots= Table[theroots[[i,2]],(i,Length[theroots]}]; 
theroots=DeleteCases[theroots, Q_Complex /; Abs[Im[Q]]>.1]; 
theroots=ReplaceAll[theroots, Q_ -> Re[Q] ]; 
rmtranscend; 

Q“.; 

theroots=Table[FindRoot[funcl==func2, (Q,theroots[[m]]} ], 

{m, Length[theroots]}]; 
theroots= Q /. theroots; 
theroots=N[Union[theroots, {}]]; 

theroots=Union[theroots, theroots + N[Table[2Pi, 
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{j, Length[theroots])]]]; 

testroots=>Union [ Select [N [ Round [ theroots* 100 ]/100 ], 

(#>=N[-Pi] 6& #<»N[(5/2) Pil)4],{)l; 
theroots=Select[theroots,(#<» N[upperlimit] kk 
#>= N[lowerlimit])6]; 

theroots=Union[N[(Round[100*theroots])/100],{)]; 

)]; 

theroots^Transpose[Select[Table[{theroots[[i]], 
(Abs[funcl-func2] /. Q-> theroots[[i]])<•01), 

{i, 1, Length[ theroots]}]/(#[[2 ] ]~=True) &]][[!]]»* 

If [theroots==*{} [ [ 1 ] ], theroots=»{ > ]; 

max={0,0}; 

For[rr=l, rr<=2, {rval=rr; ++rr; 

inax[[rval]]=Max[DeleteCases[Table[bothfunc[[rval]] 

/. Q->values, (values, 0, 2Pi, .25}], x_Complex]] }]; 

max=Max[max]; 

plotroutine; 

If[liinitsknown===n, { 

If[Length[theroots]==1, 

Print [" 

The two polar curves intersect at Q = 
topi[theroots[[1]]] ]]; 

If[Length[theroots] >1, 

Print["The two polar curves intersect at Q: 

", Map[topi,theroots] ]]; 

Pause[1];Pause[4]; 

newcurve=Input["The polar curves are plotted from Q=0 
to Q=2Pi. If you want to change the limits type y, otherwise 
type n."]; 

I f [ newcurve=”y, { 

lowerlimit=Input["Input new lower limit."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi“>Pi; 

While[lowerlimit <-N[Pi], (lowerlimit=Input["The lower limit 
can not be less than -Pi, input lower limit again."]; 
lowerlimit=lowerlimit /. sqrt->Sqrt; 
lowerlimit=lowerlimit /. pi->Pi; 

}]; 

upperlimit=Input["Input new upper limit."]; 
upperlimit=upper‘’imit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

While[upperlimit>2 N[Pi], {upperlimit=Input["The upperlimit 
can not be greater than 2Pi, input upper limit again."]; 
upperlimit=upperlimit /. sqrt->Sqrt; 
upperlimit=upperlimit /. pi->Pi; 

)]; 

rmtranscend; plotroutine 

)]; 

)]; 

]; 
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